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ABSTRACT 


Smoke-wire  flow  visualization  and  hot-wire  anemometry  have  been  used  to 


study  near  and  far  wakes  of  two-dimensional  bluff  bodies.  For  the  case  of  a  cir¬ 
cular  cylinder  at  70  <  Re  <  2000,  a  very  rapid  (exponential)  decay  of  velocity 
fluctuations  at  the  Karman  vortex  street  frequency  is  observed.  Beyond  this 
region  of  decay,  larger-scale  (lower  wave-number)  structure  can  be  seen.  In  the 
far  wake  (beyond  one  hundred  diameters)  a  broad  band  of  frequencies  is  selec¬ 
tively  amplified  and  then  damped,  the  center  of  the  band  shifting  to  lower  fre¬ 
quencies  as  downstream  distance  is  increased.  ^ 


'The  far-wake  structure  does  not  depend  directly  on  the  scale  or  frequency  of 
the  original  Karman  vortices;  the  growth  of  this  structure  is  due  to  hydro¬ 
dynamic  instability  of  the  developing  mean  wake  profile;  it  is  not  caused  by 
amalgamation  of  the  Karman  vortices.  Under  certain  conditions  amalgamation 
can  take  place,  but  is  purely  incidental,  and  is  not  the  driving  mechanism 
responsible  for  the  growth  of  larger^fcale  structure.  Similar  large  structure  is 
observed  downstream  of  porous'  flat  plates  (Re  «  6000).  which  do  not  initially 
shed  Karman-type  vortices  into  the  wake. 

^Hot-wire  measurements  show  that  two-dimensional  locally-parallel  in  viscid 
linear  stability  theory  is  adequate  to  explain  the  growth  of  downstream  struc¬ 


ture.  Namely,  measured  prominent  frequencies  In  the  cylinder  wake  are  in  close 
agreement  with  those  predicted  by  the  theory,  when  streamwise  growth  of  wake 


width  is  taken  into  account. 


Finally,  three-dimensionality  in  the  far  wake  of  a  circular  cylinder  is  briefly 


discussed. 
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Chapter  1 

Introduction 

1.1.  Historical  Overview 

Wakes  generated  by  two-dimensional  bluS  bodies  have  been  the  focus  of  hun¬ 
dreds  of  experiments  since  the  early  1900's.  Tor  a  detailed  review  or  bibliogra¬ 
phy  of  work  on  this  subject,  the  reader  is  referred  to  Morkovln  (1964)  or  Berger 
and  Wille  (1972).  The  majority  of  experimented  efforts  in  the  past  have  dealt  pri¬ 
marily  with  near  wakes  (only  a  few  diameters  downstream).  A  more  limited 
amount  of  work  has  been  done  on  the  far  wakes  of  bluff  bodies  (up  to  several 
hundred  diameters),  and  that  is  the  subject  on  which  we  wish  to  focus  here. 

Let  us  consider  the  two-dimensional  problem.  It  is  well  known  from  dimen- 

sional  analysis  and  simple  scaling  laws  that  a  plane  far  wake  grows  as  **,  while 

_  i 

the  velocity  defect  decays  as  x  T.  (See  for  example,  Cantwell  (1979)).  We  would 
expect  therefore,  as  we  go  downstream  in  the  wake,  that  the  size  of  the  Iargest- 
scale  structure  should  increase,  while  its  passage  frequency  should  decrease. 
This  in  fact  has  been  experimentally  verified.  Taneda  (1959),  using  flow  visuali¬ 
zation  in  the  wakes  of  two-dimensional  cylinders  and  fiat  plates,  observed  a  far 
wake  structure  much  like  the  original  Karman  vortex  street,  but  of  larger  scale. 
He  viewed  the  process  as  one  of  changing  hydrodynamic  stability,  l.e.,  the  origi¬ 
nal  street  decays  ("breaks  down”)  and  the  wake  "rearranges”  itself  into  a  new 
configuration  appropriate  to  the  new  (downstream)  position.  Taneda  found  that 
for  laminar  wakes  (50  <  Re  <  150),  the  ratio  of  the  secondary  vortex  street 
scale  to  that  of  the  primary  (Karman)  street  (aB/tii)  varies  from  1.5  to  3.5;  for 
turbulent  wakes  this  ratio  is  roughly  10,  and  the  secondary  street  does  not 
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appear  until  much  further  downstream,  after  the  wake  has  had  time  to  grow 
significantly.  He  also  reported  that  this  sequence  sometimes  repeats  Itself. 

Other  investigators  have  also  reported  large  scale  organized  structure  in  the 
far  wakes  of  bluff  bodies.  Grant  (195B)  observed  quasi-periodic  large  structures 
several  hundred  diameters  downstream  of  his  cylinder  at  Re  =  1300,  which  he 
described  as  “jets  of  turbulent  fluid  proceeding  outward  from  the  central  plane 
of  the  wake."  Successive  “jets"  arise  and  decay,  their  overall  size  becoming 
larger  with  increasing  downstream  distance. 

Gerrard  (1066)  studied  the  downstream  development  of  a  cylinder  wake  at  Re 
=  113.  At  x/d  =450  he  found  that  the  fundamental  (Xarman)  frequency  is  lost  in 
the  background  noise,  but  a  lower-frequency  fluctuation  is  found.  The  possibil¬ 
ity  of  a  secondary  vortex  street  was  suggested,  the  ratio  of  secondary  to  primary 
frequencies  being  about  0.033. 

Zdravkovlch  (1086)  studied  the  development  of  the  wake  behind  a  set  of  three 
circular  cylinders  in  close  proximity.  He  found  that  the  shed  vortices  rapidly 
decay  as  they  Interact  with  each  other.  Subsequently,  a  new  vortex  street  of 
larger  scale  appears  downstream.  Zdravkovlch  explained  the  formation  of  the 
secondary  vortex  street  as  a  coupling  of  two  mechanisms  -  shear  layer  instabil¬ 
ity,  and  a  rolling  up  process  Induced  by  the  distributed  vorticity  in  the  wake. 
Similar  experiments  were  done  more  recently  by  Williamson  (1963). 

An  obvious  extension  of  Zdravkovich's  work  is  the  wake  of  a  cluster  of  many 
cylinders  -  or,  in  the  general  case,  the  wake  of  a  porous  two-dimensional  body, 
such  as  a  screen  or  perforated  flat  plate  aligned  normal  to  the  flow.  Investiga¬ 
tions  of  the  latter  have  been  reported  by  Castro  (1971),  and  more  recently  by 
Valensi(l974).  Castro  categorized  two  distinct  regimes  of  flow  based  on  the  soli¬ 
dity  of  the  plate:  a)  For  high  solidity  (low  porosity)  a  Karman  vortex  street 
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dominates  the  near  wake,  b)  For  low  solidity  (c<  0.8)  the  plate  does  not  shed 
Karman  vortices,  but  there  still  Is  a  dominant  frequency  present  downstream, 
which  he  attributed  to  far  wake  instability.  Valensi  examined  the  latter  regime 
in  more  detail  with  a  53%  solid  plate;  he  shows  power  spectra  and  smoke  pic¬ 
tures.  The  dominant  near-wake  frequency  corresponds  to  a  rolling  up  of  the 
shear  layers  on  either  side  of  the  wake;  further  downstream  (beyond  five  plate 
widths),  the  dominant  frequency  Is  associated  with  a  Karman-type  vortex  street, 
which  forms  after  the  shear  layers  have  merged.  Smoke  visualization  confirms 
his  hot-wire  results. 

Another  interesting  set  of  experiments  was  reported  by  Durgin  and  Karlsson 
(1971).  They  subjected  the  cylinder  wake  to  a  deceleration,  thereby  distorting 
the  Karman  vortex  street.  Their  results  show  an  annihilation  of  the  concen¬ 
trated  vortices,  and  the  subsequent  creation  of  a  new  vortex  street  of  lower  fre¬ 
quency  and  longer  wavelength.  They,  like  Taneda,  attributed  this  effect  to 
changing  hydrodynamic  stability.  Although  the  non-decelerated  wake  of  Taneda 
does  not  experience  bending  and  stretching  of  vortex  filaments,  as  does  the 
decelerated  wake,  the  similarity  of  the  two  results  is  remarkable. 

In  a  closely  related  experiment,  Keffer  (1685)  used  a  constant-area  distorting 
duct  to  impose  a  strain  field  upon  the  wake.  His  flow  visualization  results  show 
a  secondary  vortex-street-like  structure  emerging  with  a  frequency  0.36  times 
the  Karman  frequency  at  Re  =  350. 

Interest  in  Taneda's  discovery  has  re-surfaced  in  the  past  decade.  Based  on 
an  lnvlscid  model.  Weihs  (1973)  suggested  that  multiple  Karman  vortex  street 
modes  are  possible.  In  an  attempt  to  compare  his  model  with  experiment,  he 
plotted  Taneda's  results  and  obtained  the  fit  (for  60  <  Re  <  ISO): 
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although  there  is  considerable  scatter  in  the  data.  Matsui  and  Okude  (1980) 
made  careful  hot-wire  measurements  at  Reynolds  numbers  between  100  and 
160,  of  the  velocity  fluctuations  associated  with  this  secondaLry-street 
phenomenon.  They  also  reported  an  inverse  Reynolds  number  relation,  but 
their  constant  was  much  different  than  Taneda's  (Eq.  1.1): 


nt  295 
ax  ~  Re 


(1.2) 


Thus  it  seems  that  the  development  of  the  secondary  vortex  street  is  not 
Independent  of  experimental  facility.  (Taneda’s  cylinders  were  towed  in  a  water 
channel,  while  Matsui  and  Okude  used  a  low-speed  wind  tunnel). 

Meanwhile  Townsend  (1979).  using  multiple  hot-wire  arrays,  found  periodic 
flow  patterns  resembling  Karman  vortex  streets  for  a  much  higher  Reynolds 
number  (Re  =  8000  at  xM  =  170).  "Groups”  of  three  to  five  vortical  structures 
were  observed;  the  passage  frequency  within  a  group  being  uniform,  but  varying 
considerably  from  one  group  to  the  next.  Regions  of  random  fluctuations  seem 
to  All  the  gaps  between  successive  groups.  Long-time-averaged  spectra  there¬ 
fore  do  not  show  a  peak  at  any  particular  frequency.  Townsend  estimated  that 
these  eddy  groups  contribute  15  to  20%  of  the  total  turbulent  energy.  Such 
grouping  of  structures  had  previously  been  suggested  by  Gupta  et  al.  (1971)  for 
the  turbulent  boundary  layer. 

Matsui  and  Okude  (1981)  later  abandoned  Taneda’s  original  conjecture  that 
the  secondary  street  results  from  hydrodynamic  instability,  following  the 


"breakdown”  of  the  primary  street.  In  its  stead  they  proposed  that  pairing  is 
the  mechanism  for  the  change  in  scale.  The  possibility  of  amalgamation  of 
smaller  structures,  as  in  mixing  layers,  had  also  been  speculated  by  Roshko 
(1978).  Matsui  and  Okude’s  conclusion,  however,  was  based  partly  on  flow  visual¬ 
ization.  where  the  flow  tracer  was  introduced  upstream  of  the  cylinder.  As  will 
be  pointed  out  in  Chapter  3,  this  type  of  visualization  can  be  very  misleading. 
Furthermore,  with  a  pairing  mechanism  the  frequency  of  the  secondary  street 
should  be  half  that  of  the  primary  street,  but  this  is  not  the  case  for  all  Rey¬ 
nolds  numbers.  Matsui  and  Okude  (1980)  measured  the  frequencies  of  the  pri¬ 
mary  and  secondary  streets.  From  their  Figure  10  one  can  calculate  the  ratio  of 
these  two  frequencies:  ft/fi  ranges  from  0.2  to  0.52.  Interestingly,  at  Re  =  140 
this  ratio  is  approximately  0.44;  yet  this  was  the  Reynolds  number  chosen  to 
support  their  pairing  hypothesis  (Matsui  and  Okude  (1SB1)). 

At  about  the  same  time,  Cimbala,  Naglb  and  Roshko  (19B1)  reported  strong 
experimental  evidence  in  support  of  Taneda’s  stability  hypothesis.  Namely,  they 
employed  the  smoke-wire  flow-visualization  technique  of  Corke,  et  al.  (1977)  and 
hot-wire  anemometry  to  show  that  concentrated  vortices  of  the  Karman  vortex 
street  decay  exponentially.  As  will  be  discussed  in  Chapter  3.  the  very  rapid 
decay  of  the  Karman  vortex  street  has  been  known  for  quite  some  time  (see  for 
example  Roshko  (1953)):  the  full  impact  of  this  decay,  however,  has  not  been 
realized  due  to  misinterpretation  of  streakline  flow  visualization.  Cimbala,  et  al. 
showed  that  by  100  to  150  diameters  downstream  of  a  cylinder  at  Re  »  100  ,  the 
initial  vortex  street  can  no  longer  be  found.  Thus,  they  concluded  that  the 
development  of  any  subsequent  structure  can  not  be  directly  dependent  upon 
the  original  vortex  street  (i.e.,  amalgamation),  but  rather  must  be  the  result  of 
a  local  wake  instability  of  the  developing  mean  velocity  profile.  Details  of  these 
and  subsequent  experiments  are  reported  in  the  present  paper.  Recall  that  a 
similar  conclusion  had  been  reached  by  Durgin  and  Karlsson  (1971)  for  the  case 
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of  a  decelerated  wake. 

Matsui  and  Okude  (1083),  again  studied  the  wake  at  Re  =  140.  Using  simul¬ 
taneous  hot-wire  anemometry  and  smoke  visualization  at  ar/tf=75,  they  show  evi¬ 
dence  that  the  smoke  lumps  are  indeed  vortices  at  this  location;  therefore  they 
conclude,  based  on  flow  visualization  of  these  lumps,  that  pairing  is  the  predom¬ 
inant  mechanism  for  the  change  of  scale.  The  frequency  ratio  f  z/f  \  is  not 
always  exactly  one-half  because  not  all  of  the  vortices  pair  up;  a  few  are  left 
stranded,  causing  some  irregularity.  However,  when  Matsui  and  Okude  acousti¬ 
cally  excite  the  wake  at  or  —  of  the  Karman  frequency,  they  observe  regular 
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pairing  or  tripling,  respectively. 

In  contrast  to  the  conclusions  of  Matsui  and  Okude,  Nagib  and  Desruelle 
(1983)  and  later  Desruelle  (1983)  confirmed  the  findings  of  Cimbala.  et  al.  (1981) 
with  a  different  flow  facility.  In  addition  they  used  acoustic  excitation  to  experi¬ 
mentally  derive  the  stability  characteristics  of  the  wake.  Namely,  amplification 
rates  are  plotted  as  a  function  of  disturbance  frequency  and  of  downstream 
location.  As  expected,  the  band  of  amplified  frequencies  shifts  to  lower  values 
and  broadens  as  downstream  distance  is  increased. 

Several  theoretical  and  numerical  investigations  of  this  matter  have  also 
been  made;  for  example,  see  Weihs  (1973),  Aref  and  Siggia  (1981),  Schatzman 
(1981),  and  Saffman  and  Schatzman  (1981).  For  the  most  part  these  models  are 
invlscid,  two-dimensional,  and  assume  an  infinite  array  of  vortices;  only  the 
subharmonic  instability  has  been  given  much  attention,  although  Saffman  and 
Schatzman  discuss  other  possibilities.  Robinson  (1984)  has  investigated  three- 
dimensional  stability  of  an  initially  two-dimensional  row  of  staggered  vortices. 


1.2.  Objective 

For  the  case  of  a  plane  free  shear  layer,  amalgamation  of  small  vortices  into 
larger  ones  is  an  important  contributor  to  the  growth  of  coherent  structures. 
(See  for  example  Brown  and  Roshko  (1674).)  To  what  extent  amalgamation 
influences  the  scale  of  structure  in  a  luaJea  is  of  concern  here.  As  discussed 
above,  there  Is  considerable  disagreement  among  researchers  as  to  the 
mechanism(s)  responsible  for  the  large-scale  structure  observed  in  far  wakes  of 
two-dimensional  bluff  bodies.  The  two  views  are  summarized  here: 

a)  Hydrodynamic  instability  of  the  mean  wake  profile  is  responsible  for  the 
growth  of  downstream  structure,  independent  of  and  not  directly  resulting  from 
the  vortices  shed  from  the  body  (Karman  vortex  street). 

b)  The  large  structure  in  the  far  wake  is  a  result  of  amalgamation  of  the  Kar¬ 
man  vortices  into  larger  vortical  structures. 

Our  overall  objective  then  is  to  describe  definitively  the  mechanism  that 
causes  the  change  of  scale.  The  experimental  tools  we  have  employed  toward 
this  end  are  hot-wire  anemometry  and  smoke-wire  visualization.  We  are  pri¬ 
marily  concerned  with  the  two-dimensional  aspect  of  the  problem:  three- 
dimensionality  is  also  discussed  briefly. 


Chapter  2 


Experimental  Facility  and  Instrumentation 

2.1.  Facility 

The  flow-visualization  experiments,  as  well  as  the  hot-wire  measurements, 
were  conducted  in  an  open-return,  low-turbulence  wind  tunnel,  shown  schemati¬ 
cally  in  Figure  2.1.  The  test  section  was  0.5  m  square  and  2  m  long;  flow  velocity 
was  variable  from  0.5  m/  s  to  12  m/  s.  The  turbulence  intensity  u'rma/Um  of  the 
fre68tream  was  originally  measured  to  be  0.4%  at  Um  =  2  m/  s.  We  had  con¬ 
cluded  that  this  turbulence  level  was  too  high  for  accurate  wake  measurements, 
where  the  fluctuations  in  the  wake  itself  are  quite  low.  Therefore  we  refurbished 
the  facility  with  a  new  test-section  and  added  several  screens  in  the  settling 
chamber,  in  the  manner  suggested  by  Loehrke  and  Nagib  (1972).  The  final 
result  in  the  test  section  is  a  background  turbulence  level  of  about  0.1  to  0.2%  at 
Um  -  2  m/  s. 

2.2.  Hodeks 

The  circular  cylinders  used  in  the  experiments  varied  in  diameter  from  0.06 
cm  to  0.5  cm,  giving  a  Reynolds  number  range  from  70  to  over  2000  (based  on 
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cylinder  diameter.  Re  =  — — — ).  The  cylinders,  made  of  drill  rod,  spanned  the 

test  section  and  passed  through  the  walls  at  either  end.  where  they  were  held  in 
tension  by  fasteners  outside  the  tunnel  walls.  Drill  rod  provided  the  required 
strength,  uniformity  of  diameter,  and  smoothness.  Before  each  run,  the 
cylinder  was  polished  carefully  to  remove  any  dust  build-up  or  corrosion, 

Three-dimensionally  perturbed  cylinders  were  used  as  well  as  the  smooth  two- 
dimensional  ones,  (see  Chapter  6).  To  make  one  of  these  three-dimensional 


cylinders,  paint  was  sprayed  on  the  surface  of  the  cylinder  at  specific  intervals, 
and  for  certain  lengths  along  its  axis,  after  masking  those  sections  which  were 
to  remain  unpainted.  In  this  manner,  the  cylinder  diameter  was  increased 
where  painted,  producing  a  three-dimensional  disturbance  along  the  span. 

In  addition  to  circular  cylinders,  porous  flat  plates  were  also  stretched  across 
the  test  section,  with  their  broad-sides  facing  the  flow.  The  plates  were  approxi¬ 
mately  2  cm  in  width  and  were  cut  from  screens.  (Reynolds  number  based  on 
this  width  ranged  from  3000  to  9000.)  In  the  initial  stages  of  this  work,  done  in 
cooperation  with  H.  Nagib,  we  looked  at  a  number  of  porous  plates  of  various 
solidities  and  mesh  sizes.  From  these  we  picked  one  with  initially  laminar  flow, 
and  one  with  initially  turbulent  flow  for  our  smoke-visualization  experiments 
(see  Chapter  4):  a)  a  29%  solidity  screen  with  wire  diameter  0.18  mm  and  mesh 
size  1.15  mm.  and  b)  a  47%  solidity  screen  with  wire  diameter  0.43  mm  and 
mesh  size  1.59  mm,  respectively.  The  screens  were  held  in  tension  by  fasteners 
which  passed  through  the  tunnel  walls. 

2.3.  Smoke-Wire  Flow-Visualization  Technique 

A  very  simple  technique  for  introducing  closely  spaced  streaklines  in  wind 
tunnels  is  called  “smoke-wire  flow  visualization.”  and  is  described  by  Corke,  et 
al.  (1977).  In  the  present  experiments,  the  smoke-wire  used  for  flow  visualiza¬ 
tion  was  a  0.13  mm  diameter  stainless  steel  wire  which  stretched  from  top  to 
bottom  of  the  test  section  through  the  tunnel  walls.  The  top  end  of  the  wire 
passed  through  a  No.  27G  hypodermic  needle  and  then  into  a  3.2  mm  hollow 
brass  cylinder,  which  was  clamped  to  the  outside  of  the  tunnel.  This  cylinder 
served  a  two-fold  purpose:  a)  to  hold  the  smoke-wire  in  place,  and  b)  to  serve  as 
a  channel  for  the  oil  to  reach  the  wire,  as  described  below.  The  other  end  of  the 
smoke-wire  passed  through  the  bottom  of  the  test  section  where  it  was  clamped 
to  a  hanging  weight  which  provided  constant  tension  in  the  wire.  At  either  end 


electrical  connections  were  made  so  that  current  could  pass  through  the  wire. 
The  scenario  is  as  follows:  Mineral  oil  is  pumped  through  the  hollow  cylinder  and 
through  the  hypodermic  needle,  forcing  a  large  drop  to  fall  down  the  smoke- 
wire.  As  the  drop  falls,  surface  tension  causes  the  oil  along  the  wire  to  break  up 
into  tiny  beads  or  droplets.  Two  or  three  initial  drops  are  typically  used  to  pro¬ 
vide  a  uniform  distribution  of  smaller  droplets  along  the  span  of  the  smoke- 
wire.  At  this  point  a  current  is  passed  through  the  wire,  which  heats  it  up 
sufficiently  to  cause  the  oil  to  burn.  The  result  is  rows  of  streaklines  which 
mark  the  flow  downstream  of  the  smoke-wire.  (The  amount  of  current  needed 
to  accomplish  this  varies  with  tunnel  speed,  but  is  typically  one  ampere  at  a 
speed  of  2  m/  s.)  After  an  adjustable  time  delay  to  allow  the  smoke  to  convect 
completely  into  the  field  of  view,  a  camera  is  electronically  triggered  along  with 
a  strobe  flash  to  record  the  streakline  pattern  on  film. 

The  electronic  synchronizing  controller  for  the  smoke-wire  operation  was 
built  at  the  Illinois  Institute  of  Technology,  where  this  technique  was  perfected. 
For  further  details  of  its  usage  see  Corke,  et  al.  (1977).  A  Pentax  ME-Super  35 
mm  camera  with  motor  drive  was  used,  along  with  a  General  Radio  Model  1540 
Stroboscope.  We  had  the  best  success  with  Kodak  Tri-x  film  pushed  one  stop, 
and  printed  on  F-5  high  contrast  paper.  The  test  section  of  the  tunnel  was  plexi¬ 
glass  on  3  sides  (top,  bottom  and  front),  and  the  back  was  wood  that  had  been 
spray-painted  a  very  flat  black.  The  strobe  was  aligned  perpendicular  to  the 
flow,  shining  up  through  the  bottom,  and  the  camera  was  positioned  about  a 
meter  from  the  front  of  the  test  section. 

2.4.  Instrumentation 

Freestream  velocities  were  measured  with  a  United  Sensors  pitot-static  tube 
and  am  electronic  Barocel  manometer.  The  pitot-static  tube  was  placed  half-way 
between  the  top  of  the  test  section  and  the  model  axis  in  order  to  obtain 


accurate  ireestream  measurements  as  well  as  to  minimize  interference  with  the 
flow  over  the  model. 

Fluctuating  velocities  were  measured  with  miniature  hot-wire  probes,  pur¬ 
chased  from  T.S.l.  (models  1260-T1.5  and  1248-T1.5).  The  single-sensor  probe 
was  a  4  /xm  platinum-plated  tungsten  wire,  with  a  length  of  1.5  mm  {L/d  =  375). 
It  was  mounted  parallel  to  the  longitudinal  axis  of  the  model,  and  with  the  probe 
body  oriented  at  45°  with  respect  to  the  x  -axis  (direction  of  freestream).  A  six 
mm  cylindrical  strut  supported  the  hot-wire  assembly  from  the  top;  it  was 
covered  with  an  extruded  aluminum  airfoil-shaped  tube  to  minimize  vibration 
and  disturbance.  The  entire  assembly  was  mounted  to  am-y  traverse  which 
rested  on  top  of  the  test  section.  The  dual-sensor  probe  consisted  of  two  wires 
in  nn  "x"  configuration,  each  wire  45°  with  respect  to  the  freestream,  but  per¬ 
pendicular  to  its  neighbor.  The  wires  were  of  the  same  diameter  as  that  of  the 
single-sensor  probe  described  above.  The  dual-sensor  probe  was  mounted  paral¬ 
lel  to  the  freestream  such  that  the  “x”  lay  in  the  x~y  plane.  In  this  orientation 
•u'  and  v’  fluctuations  could  be  measured  as  described  in  Appendix  B.  The 
strut,  traverse,  etc.,  were  identical  to  those  described  above  for  the  single¬ 
sensor  probe. 

A  two-channel  constant-temperature  anemometer  circuit  was  built  in-house 
by  Nosenchuck,  et  al.  (1983).  Some  modifications  to  the  circuit  were  required  to 
obtain  a  very  high  slgnal-to-noise  ratio.  Its  output  was  analyzed  in  one  of  two 
ways: 

a)  For  single-sensor  real-time  spectral  analysis  an  HP  35B2A  spectrum 
analyzer  was  used,  along  with  an  x  -y  plotter. 

b)  For  dual-sensor  time-trace  measurements,  and  for  simultaneous  hot-wire 
and  smoke-wire  measurements,  a  digital  data  acquisition  system  was  used. 
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Following  the  schematic  diagram  of  Figure  2.2,  the  anemometer  output  was  sent 
to  a  signal  conditioner  which  offset  the  DC  component,  amplified  the  signal  to 
cover  the  full  range  of  the  12-bit  A/D  converter  (0  -  10  Volts),  and  applied  a  low- 
pass  filter  as  an  anti-aliasing  precaution.  An  6-bit  microcomputer  was  used  to 
sample  and  record  the  data  digitally  onto  floppy  disk.  At  a  programmable  time 
during  sampling,  the  smoke-wire  could  be  turned  on,  as  well  as  the  camera  and 
strobe,  in  order  to  record  simultaneous  hot-wire  data  and  smoke-wire  photo¬ 
graphs.  The  signal  conditioner  and  relay  unit  were  built  by  the  author;  the 
microcomputer  was  also  built  in-house  by  Nosenchuck,  et  al.  (1983). 

The  hot-wire  data  were  subsequently  analyzed  on  GALCIT’s  DEC  PDP-11  com¬ 


puter  system. 
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Chapter  3 

The  Art  of  Flow  Visualisation 

3.1.  Integration  Effect  of  Streaklines 

Flow  visualization,  though  a  very  useful  tool  in  fluid  mechanics  research,  can 
often  be  misleading,  as  has  been  pointed  out  by  Hooker  (1936)  and  more 
recently  by  Hama  (1962).  The  main  problem  with  unsteady  flows  is  that  streak- 
lines  are  not  equivalent  to  streamlines,  even  though  they  are  sometimes  con¬ 
fused  as  such.  In  the  laboratory,  it  is  often  more  convenient  to  manufacture 
streaklines  than  streamlines;  hence  the  former  are  much  more  popular  among 
experimentalists.  The  most  common  method  of  streakline  visualization  is  the 
steady  introduction  of  a  flow  tracer  (typically  dye  or  smoke)  from  a  fixed  loca¬ 
tion  in  the  flow.  For  wakes  this  location  is  usually  upstream  of  the  body,  or  else 
on  the  surface  of  the  body  Itself.  But  a  dye  or  smoke  filament  is  distorted  as  it 
travels  downstream,  and  the  streakline  pattern  seen  at  some  downstream  loca¬ 
tion  contains  information  integrated  all  the  way  back  to  its  point  of  introduc¬ 
tion.  Thus,  when  studying  a  streakline  one  cannot  be  sure  if  a)  the  flow  at  a 
given  location  is  currently  experiencing  the  distortions  implied  by  the  streakline 
pattern,  or  b)  that  pattern  is  merely  a  remnant  of  the  streakline’s  past  history 
(an  integrated  "memory”  as  it  were).  Therefore,  in  order  to  visualize  the  true 
nature  of  the  flow  at  a  given  location,  it  Is  desirable  to  introduce  the  flow  tracer 
as  close  to  that  location  as  possible.  Traditionally  the  hydrogen  bubble  tech¬ 
nique  in  water  has  been  useful  because  of  this  advantage;  the  smoke-wire  tech¬ 
nique  in  air  is  analogous  and  has  been  employed  here. 

As  the  case  in  point  consider  the  two-dimensional  Karman  vortex  street  in 
the  wake  of  a  circular  cylinder  at  low  Reynolds  number.  Figure  3.1  shows  a 
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series  of  photographs  taken  at  different  times  and  at  progressively  Increasing 
downstream  locations  in  the  wake  at  Re  =  90.  The  flow  is  from  left  to  right,  with 
the  cylinder  axis  perpendicular  to  the  plane  of  view.  In  each  case  the  smoke- 
wire  has  been  placed  at  the  left-most  edge  of  the  field  of  view,  and  the  photo¬ 
graphs  are  aligned  so  that  any  vertical  line  passes  through  the  same  down¬ 
stream  location  in  all  of  the  photographs.  In  photo  a)  the  smoke-wire  is  just 
downstream  of  the  cylinder;  the  vortex  street  is  clearly  marked  and  remains 
visible  to  the  downstream  end  of  the  photograph  ( z/d  -  200).  From  photo  a) 
alone  one  might  conclude  that  the  Karman  vortex  street  is  still  active  at  this 
downstream  location.  In  photos  b)  and  c)  the  smoke-wire  has  been  moved  down¬ 
stream  to  50  and  100  diameters,  respectively.  The  vortices  are  Been  to  be  much 
weaker  for  these  smoke-wire  positions  than  would  appear  from  photograph  a) 
alone.  In  photo  d)  the  smoke-wire  is  at  xM  -  150  and  all  the  streaklines  are 
parallel,  which  implies  that  the  street  at  this  location  is  so  weak  it  can  not  be 
detected  by  the  smoke-wire  technique.  This  essentially  parallel  wake  profile  con¬ 
tinues  downstream  with  no  further  evidence  of  discrete  vortices  all  the  way  to 
350  diameters. 

The  integration  effect  of  streaklines,  as  described  above,  is  clearly  demon¬ 
strated  by  comparing  the  flow  pattern  observed  at  z/tf  =  200  in  each  of  the  pho¬ 
tographs  in  Figure  3.1;  the  pattern  is  drastically  different,  depending  on  smoke- 
wire  position.  Photo  a)  shows  what  appears  to  be  the  Karman  vortex  street  at 
z/tf  =  200.  Photos  b)  and  c)  show  a  similar  structure  but  much  weaker,  while 
photo  d)  indicates  that  no  structure  is  present  at  this  wake  location.  Which 
photo  then  can  we  trust?  That  is,  how  can  we  objectively  describe  the  flow  field 
from  this  confusing  array  of  photographs?  It  is  our  contention  that  in  order  to 
accurately  discern  the  flow  at  some  location,  the  smoke-wire  must  be  placed  at  a 
proper  distance  upstream  of  that  location.  If  it  is  placed  too  far  upstream,  the 
streaklines  become  very  distorted  and  may  lead  us  to  erroneous  conclusions. 
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likewise  if  it  is  placed  too  close  to  our  observation  point,  the  streaklines  may 
not  have  had  time  to  deform,  again  leaving  us  confused. 

The  clue  to  understanding  this  figure  lies  in  the  following:  Suppose  a  streak¬ 
line  rolls  up  as  it  marks  a  vortex  moving  downstream.  The  resulting  pattern  will 
clearly  indicate  that  vortex.  Now  suppose  that  the  vortex,  as  it  travels,  decays 
at  a  much  faster  rate  than  the  smoke  diffuses.  We  can  then  argue  that  the 
smoke  remains  in  its  integrated  pattern,  even  after  the  vorticity  has  largely 
diffused  (leaving  the  vortex  extremely  weak).  The  streakline  pattern  is  now 
fixed,  and  is  simply  convectlng  along  with  the  mean  flow,  no  longer  changing  its 
shape.  This  is  exactly  what  we  observe  in  photo  a)  of  Figure  3.1;  notice  that 
beyond  100  diameters  or  so,  the  steakline  pattern  no  longer  changes 
significantly  —  it  merely  convects  downstream.  There  is  of  course  some  shearing 
of  the  pattern,  as  the  outside  of  the  wake  convects  faster  than  the  center.  This 
can  be  seen  in  photo  b)  as  well  as  in  photo  a)  of  Figure  3.1.  Zdravkovich  (1969) 
shows  what  he  describes  as  extremely  elongated  elliptical  vortices  at  x/tf  »  600; 
his  smoke  was  introduced  at  the  cylinder,  and  therefore  suffers  the  same 
integration  effect  as  our  photo  a). 

Keeping  these  concepts  in  mind,  we  can  now  answer  the  questions  posed 
above.  Each  photograph  in  the  figure  is  useful  only  for  a  finite  distance  down¬ 
stream  of  the  smoke-wire.  To  obtain  an  accurate  description  of  the  entire 
developing  flow  field,  it  is  necessary  to  place  the  smoke-wire  at  various  positions, 
and  then  look  at  the  collage  of  pictures  as  a  whole.  For  example,  the  very  rapid 
decay  of  the  Karman  vortex  street  is  not  clear  from  any  one  of  the  photographs 
in  Figure  3.1;  but  when  the  entire  figure  is  studied  simultaneously,  the  decay  is 
quite  evident.  We  should  mention  here  that  the  smoke-wire  itself  has  a  wake, 
but  the  Reynolds  number  based  on  wire  diameter  Is  only  about  seven.  At  Rey¬ 
nolds  numbers  below  forty  or  so,  the  wake  of  the  smoke-wire  is  steady  and  only 


-  la  - 


influences  the  flow  immediately  downstream.  We  therefore  consider  the  smoke- 
wire  as  essentially  non-intrusive. 

At  this  point  the  reader  may  be  questioning  the  validity  of  the  above  argu¬ 
ment  for  air  flow,  where  the  ratio  of  viscous  to  molecular  diffusivity  (Schmidt 
number)  is  of  order  unity.  One  may  argue  that  as  a  structure  in  the  flow  is  dis¬ 
sipated  by  viscosity,  the  smoke  particles  are  diffused  at  nearly  the  same  rate; 
smoke  should  not  remain  in  a  pattern  for  any  longer  time  period  than  the  life¬ 
time  of  the  structure  generating  that  pattern.  A  simple  explanation  concerning 
this  apparent  dilemma  is  as  follows:  Smoke  in  air  is  actually  composed  of  tiny 
aerosol-type  particles,  whose  mass  is  significantly  greater  than  the  mass  of  sur¬ 
rounding  air  molecules.  Hence  diffusion  of  these  particles  is  extremely  slow  on  a 
molecular  scale;  we  may  then  define  an  “effective  Schmidt  number"  which  is 
orders  of  magnitude  larger  than  that  for  air  alone.  Lapple  (1961)  has  calculated 
the  diffusion  coefficient  of  particulate  matter  into  air  at  25"  C.  For  oil  smokes, 
particle  size  is  approximately  0.2  /xm,  which  gives  a  diffusion  coefficient  of  about 
2  x  10"®  cmz/s.  In  air,  v  *  0.l6cm*/s;  Our  effective  Schmidt  number  is  thus  of 
order  10s. 

That  the  diffusion  of  smoke  is  indeed  slow  is  clearly  seen  in  Figure  3.1;  if  one 
follows  a  single  streakline  in  the  freestream  outside  of  the  wake,  no  significant 
dispersion  of  that  streakline  is  observed  (i.e.,  it  does  not  broaden  as  it  travels 
downstream).  We  can  quantitatively  estimate  (by  dimensional  reasoning)  that 
the  width  of  a  streakline  in  the  freestream  should  grow  like 
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where  D  is  the  mass  diffusivity  and  w  the  width  of  the  streakiine.  In  Figure 
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3.1,  Um  =  1 15  cm/  s,  d  =  0.13  cm,  and  the  total  distance  in  our  field  of  view  (200 
diameters)  corresponds  to  about  25  cm.  For  D  *  2  x  10®  cm'/s,  the  expected 
width  increase  at  the  end  of  the  photograph  is  Aw  *>  6.6  x  10-4  cm  (5  x  10~® 
diameters).  Such  a  small  width  increase  would  be  hardly  noticeable;  indeed  in 
Figure  3.1  the  initial  width  of  the  streakline  is  smaller  than  one  cylinder  diame¬ 
ter.  and  does  not  change  significantly  downstream. 

3.2.  Numerical  Simulation  of  Streaklines 

It  is  a  relatively  simple  exercise  to  model  streakline  patterns  numerically  for 
a  given  velocity  field.  Following  Hama  (1982),  we  introduce  marked  "particles'’ 
at  a  fixed  point,  and  track  each  particle  in  time  as  it  is  carried  away  by  the  velo¬ 
city  field  (Lagrangian  reference  frame).  Plotting  a  curve  through  each  particle 
results  in  a  numerical  streakline  (of  infinite  Schmidt  number).  This  method  was 
adapted  here  to  simulate  smoke-wire  streaklines  numerically  at  various  down¬ 
stream  positions  in  a  plane  wake. 

The  velocity  field  chosen  for  this  simulation  is  a  convecting  ideal  Karman  vor¬ 
tex  street  (Lamb  1945),  modified  to  include  the  effect  of  viscosity.  This  model 
was  originally  suggested  by  Hooker  (1936),  and  recently  has  been  used  by  Davies 
(1976);  details  are  given  in  Appendix  A. 

Figure  3.2  shows  numerical  results  which  simulate  a  cylinder  wake  at  He  =  90. 
The  “smoke-wire"  is  located  at  xM  =  0  and  at  xM  -  100.  Comparison  of  the 
two  patterns  shows  the  misleading  integration  effect  of  streaklines.  When  the 
point  of  introduction  of  the  flow  tracer  is  far  upstream  of  the  wake  location 
being  considered,  its  entire  history  is  contained  in  the  streakline  pattern,  and 
may  confuse  the  intepretation.  The  conclusion  from  this  exercise  again  is  that 
extreme  care  must  be  taken  when  interpreting  streakline  patterns! 
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Chapter  4 

Flow  Visualization  Results 

4.1.  Decay  of  the  Karman  Vortex  Street 

As  discussed  in  Chapter  3,  an  observation  immediately  apparent  from  Figure 
3.1  is  the  very  rapid  decay  of  the  Karman  vortex  street.  This  decay  was  found  to 
occur  at  all  Reynolds  numbers  we  investigated.  In  each  case,  when  the  smoke- 
wire  is  placed  far  enough  downstream  of  the  cylinder,  the  Bhed  vortices  can  no 
longer  be  found.  Furthermore,  the  location  of  their  disappearance  seems  to 
move  upstream  with  Reynolds  number.  For  Re  =  90  it  occurs  at  about  125  to 
150  diameters.  For  Re  =  155,  where  the  wake  is  just  on  the  verge  of  transition  to 
turbulent  flow,  the  Karman  vortex  street  can  not  be  visualized  beyond  75  to  100 
diameters.  (Of  course,  in  this  case  a  secondary  wake  structure  is  rapidly 
developing  at  this  downstream  location,  as  will  be  discussed  in  the  following  sec¬ 
tion.)  For  transitional  and  turbulent  wakes  (Re  >  160)  the  vortices  are  dissi¬ 
pated  much  more  rapidly,  and  are  not  seen  beyond  xM  =50. 


The  rapid  decay  of  the  Karman  vortex  street  has  been  known  for  quite  some 
time.  For  example  Roshko  (1953)  found  that  for  Re  =  50  and  100.  the  energy 
intensity  of  the  Karman  vortices  is  so  low  beyond  100  diameters,  that  back¬ 
ground  tunnel  turbulence  becomes  significant.  For  Reynolds  numbers  in  the 
“transition  range"  (150  <  Re  <  300),  the  decay  is  even  more  pronounced.  For 
turbulent  wakes  Roshko’s  results  are  again  supported  by  the  present  experi¬ 
ments,  where  Karman  vortices  could  not  be  seen  beyond  50  diameters. 
Nevertheless,  it  seems  that  the  full  implication  of  this  decay,  particularly  for 
laminar  vortex  streets,  has  been  stifled  due  to  misinterpretation  of  streakline 
photographs.  For  example,  Zdravkovich  (1989)  shows  what  appear  to  be  very 
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elongated  K&rman  vortices  at  x/d  =  600!  In  light  of  our  discussion  of  streakline 
integration  effects  in  Chapter  3,  his  elongated  ellipses  are  merely  the  integrated 
pattern  from  far  upstream.  Tritton  (1977)  states  that  "...  the  vortex  street  con¬ 
tinues  to  all  distances  downstream."  again  a  misinterpretation  of  streakline  pat¬ 
terns.  In  actuality,  the  Karman  vortex  street  has  completely  decayed  by  at  most 
one  or  two  hundred  diameters.  Indeed,  based  on  this  rapid  decay,  the  rows  of 
vortices  should  not  really  be  called  a  “street"  at  all  (the  original  term  having 
implied  a  long,  regularly-spaced  pattern).  In  Chapter  5  we  take  a  quantitative 
look  at  this  decay,  using  hot-wire  anemometry. 

4.2.  The  Secondary  Vortex  “Street" 

Because  we  are  able  to  mount  smoke-wires  at  any  desired  position,  the 
smoke-wire  technique  is  convenient  for  studying  development  of  downstream 
structure  in  a  wake.  In  particular,  placement  of  a  smoke-wire  at  positions 
beyond  which  the  Karman  vortex  street  has  largely  decayed  yields  some 
interesting  results.  For  Reynolds  numbers  less  than  about  100,  we  do  not 
observe  any  further  reorganization  or  development  of  large  structure  in  the  far 
wake.  Figure  3.1  shows  a  typical  case  (Re  =  90).  Notice  that  beyond  xM  =  150, 
the  streaklines  in  photo  d)  are  parallel,  with  no  apparent  regions  of  concen¬ 
trated  vorticity.  For  these  low-Reynolds-number  wakes,  large  secondary  struc¬ 
ture  could  not  be  seen  as  far  downstream  as  we  investigated,  which  was  about 
500  diameters  (not  shown  in  Figure  3.1). 

For  100  <  Re  <  160,  a  structure  quite  similar  to  that  of  the  Karman  vortex 
street,  but  of  larger  scale,  is  observed.  Three  cases  in  this  Reynolds  number 
range  are  shown  in  Figures  4.1  through  4.3,  at  Re  =  130,  140,  and  155  respec¬ 
tively.  The  photos  are  aligned  in  the  same  manner  described  for  Figure  3.1.  The 
“strength"  of  the  secondary  street  may  be  determined  by  how  quickly  the 
streaklines  deform.  From  these  photographs  and  others,  we  observe  that  the 
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strength  increases  with  Reynolds  number  in  the  range  100  <  Re  <  160.  The 

» 

scale  of  this  structure  is  roughly  2  to  3  times  that  of  the  primary  Karman  street. 

Figure  4.3  shows  some  of  the  strongest  and  most  regular  secondary  structure 
that  we  have  recorded  on  film.  Photos  b)  and  c)  represent  identical  tunnel  con¬ 
ditions,  but  were  taken  about  ten  seconds  apart.  Both  are  included  here  to  illus-  * 

trate  that  the  secondary  structure  is  not  always  continuous  -  it  comes  and  goes 
in  time.  Photo  b)  shows  an  irregular  secondary  street;  this  is  a  more  typical  case 
than  photo  c),  where  an  unusually  regular  and  strong  street  has  been  captured.  » 

Three-dimensionality  may  play  an  important  role  in  determining  how  regular 
the  downstream  structure  appears  in  these  edge  views.  We  discuss  three- 
dimensionality  in  Chapter  8. 

For  Re  >  160  the  wake  becomes  turbulent,  making  flow  visualization 
extremely  difficult.  Taneda  (1959)  reports  having  seen  a  secondary  street  at 
high  Reynolds  number,  but  only  after  several  hundred  diameters.  Matsui  and  j 

Okude  (1980)  could  not  And  any  such  street  for  Re  >  160,  although  they  had  not 
looked  as  far  downstream  as  had  Taneda.  In  the  present  experiments  we  did  not 
find  a  well-organized,  easily-recognizable  secondary  "street”  for  Re  >  180.  How- 

I 

ever,  "groups"  or  "bursts"  of  large  vortical  structures  can  sometimes  be 
observed,  similar  to  those  observed  by  Grant  (1958)  and  Townsend  (1979).  Fig¬ 
ure  4.4,  for  example,  shows  a  cylinder  wake  at  Re  =  190.  Notice  that  after  the 
Karman  street  decays,  the  wake  becomes  very  disorganized.  However,  if  we  follow  * 

a  streakline  in  the  outer  edge  of  the  wake  at  x/d.  **  200,  (photo  b))  we  see  a 
somewhat  regular  wavy  pattern,  of  a  scale  much  larger  than  the  original  Kar¬ 
man  street.  Figure  4.5  shows  the  wake  at  Re  =  2200;  some  of  this  structure  can  I 

also  be  observed  in  photo  b)  at  z/d  *  250.  Here  the  structure  is  more  regular 
than  for  the  case  at  Re  =  190. 
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4.3.  Porous  Plat  Plate  Takes 

It  is  important  to  note  that  in  ail  the  above  cases,  secondary  structure 
emerges  after  the  primary  (Karman)  street  has  largely  decayed.  The  secondary 
structure  therefore  appears  to  develop  independently;  for  example  its  scale  (or 
frequency)  is  not  necessarily  in  a  fixed  ratio  to  that  of  the  primary  street.  This 
observation  leads  us  to  suspect  that  the  secondary  structure  results  from  an 
instability  of  the  mean  wake  profile.  It  therefore  seems  reasonable  to  assume 
that  if  we  can  produce  a  wake  which  initially  has  no  Karman  vortex  street,  a 
street-like  structure  may  emerge  downstream  due  to  wake  instability.  Such  a 
wake  is  produced  by  a  porous  flat  plate  aligned  normal  to  the  flow  direction, 
provided  the  solidity  a  is  lower  than  about  80%,  according  to  Castro  (1971). 


Two  porous  plate  wakes  are  shown  in  Figures  4.6  and  4.7.  In  the  first  case, 
the  plate  is  very  porous  (  a  =  29%),  and  the  wake  immediately  downstream 
shows  no  structure  of  any  scale.  Presumably  all  that  exists  is  a  laminar  "top- 
hat"  velocity  profile.  At  about  10  diameters  downstream  (Figure  4.6a),  a  shear 
layer  instability  leads  to  the  growth  of  small  vortices  on  either  side  of  the  wake. 
These  vortices  amalgamate,  and  as  they  grow  they  arrange  themselves  eventu¬ 
ally  into  a  fairly  regular  vortex  street  pattern  at  x/d  «  20,  presumably  caused 
by  far-wake  instability.  At  conditions  similar  to  ours,  Castro  (1971)  observed  a 
dominant  frequency  in  the  wake,  which  he  also  attributed  to  hydrodynamic  ins¬ 
tability  of  the  mean  velocity  profile. 

The  second  case  (Figure  4.7)  is  for  a  porous  flat  plate  with  a  much  higher  soli¬ 
dity  (  a  =  47%),  but  at  approximately  the  same  Reynolds  number.  The  obvious 
difference  here  is  the  abundance  of  small-scale  turbulent  structure  immediately 
downstream  of  the  plate.  Notice  however  that  the  body  is  still  not  solid  enough 
to  shed  a  Karman  vortex  street.  Nevertheless,  this  wake  also  goes  unstable 
downstream,  and  by  30  or  40  diameters  a  well-defined  vortex  street  is  observed. 


Indeed,  if  we  could  ignore  the  small-scale  structure  superposed  on  this  flow,  Fig¬ 
ure  4.7  is  strikingly  similar  to  Figure  4.6.  Valensi  (1974)  shows  a  similar  photo¬ 
graph  for  a  plate  with  a  =  53%  at  Re  =  5  x  104. 

Hot-wire  measurements  for  the  second  case  above  are  discussed  in  Chapter  5. 
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Chapter  5 

Hot-Wire  Measurements 

5.1.  Exponential  Decay  of  the  Karmen  Vortex  Street 

Smoke-irire  flow  visualization  of  the  far  wakes  of  bluff  bodies  has  given  us  new 
qualitative  insight  into  Taneda’s  phenomenon  of  vortex  street  breakdown  and 
rearrangement.  We  now  present  the  results  of  hot-wire  measurements  which 
give  us  quantitative  information  about  this  process. 

To  obtain  detailed  velocity  measurements  in  the  wake,  we  used  a  "cross-"  or 
"x-"  wire,  as  discussed  in  Chapter  2,  along  with  a  digital  data  acquisition  system. 
The  calibration  procedure  and  details  of  the  data  acquisition  are  described  in 
Appendix  B.  Specific  details,  such  as  freestream  velocity,  cylinder  diameter, 
mean  velocity  profiles,  etc.,  as  well  as  a  description  of  the  methods  used  to 
analyze  the  data,  can  also  be  found  in  that  appendix. 

First,  let  us  consider  the  decay  of  the  ICarman  vortex  street  as  downstream 
distance  is  increased.  As  had  been  pointed  out  by  Bevilaqua  (1975),  velocity 
fluctuations  due  to  a  vortex  street  are  more  readily  detected  in  the  cross-wake 
( v )  velocity  component.  Hence  consider  the  v'  fluctuations,  where  v  =  V  +  v’  . 
In  Figure  5.1  we  plot  the  amplitude  spectrum  of  v'  /Um  for  a  circular  cylinder 
wake  at  Re  =  150  and  at  xM=  25.  (Amplitude  spectrum  is  the  square  root  of 
power  spectrum;  the  units  are  in  percent  of  freestream  velocity.)  As  expected,  a 
delta  function  spike  is  seen  at  the  shedding  frequency,  /*,  which  is  the  only 
significant  frequency  in  the  near  wake.  Such  spectra  were  taken  at  each  r- 
location,  averaged  at  the  two  y -locations  corresponding  to  inflection  points  in 
the  mean  velocity  profile  (one  on  either  side  of  the  wake).  Discrete  spectral 


amplitude  at  /  =  /*  is  plotted  against  downstream  distance  x/d  in  Figure  5.2 
for  Re  =  140  and  well  as  for  Re  =  150.  For  both  Reynolds  numbers,  the  amplitude 
is  seen  to  decay  exponentially  for  25  <  x/d  <  150.  By  least-squares  fitting  a 
straight  line  through  the  first  6  or  7  data  points  in  Figure  5.2,  we  obtain  initial 
decay  rates  of  0.0246  decades  per  diameter  for  Re  =  150,  and  0.0249  decades 
per  diameter  for  Re  =  140.  Similar  decay  rates  have  been  reported  for  the  W 
fluctuations  by  other  investigators  as  well.  Hussain  and  Ramjee  (1976)  obtained 
hot-wire  traces  of  u'  at  Re  =  145  for  various  downstream  distances.  From  their 
Figure  4,  the  amplitude  of  the  fluctuations  can  be  seen  to  decay  exponentially  in 
the  region  5  <  x/d  <  60.  In  our  earlier  work  (Cimbala,  et  al.  (1981))  at  Re  =  155, 
we  measured  discrete  u'  fluctuations  at  the  shedding  frequency,  and  found  u'^, 
to  decay  exponentially  up  to  100  diameters.  Matsui  and  Okude  (1981)  made  simi¬ 
lar  measurements  at  Re  =  140,  while  Desruelle  (1983)  showed  an  exponential 
decay  for  three  different  Reynolds  numbers.  The  decay  rates  for  all  these  cases 
are  summarized  below: 


Reference 

Velocity 

Re 

Decay  Rate 

(decade/diameter) 

Desruelle  (19B3) 

u' 

80 

0.0143 

Desruelle  (1983) 

u' 

120 

0.0169 

Present  Results 

u‘ 

140 

0.01B7 

Present  Results 

v' 

140 

0.0249 

Matsui  &  Okude  (1981) 

u’ 

140 

0.0170 

Hussain  &  Ramjee  (1976) 

u' 

145 

0.0156 

Present  Results 

u‘ 

150 

0.0209 

Present  Results 

V1 

150 

0.0246 

Desruelle  (1983) 

u' 

155 

0,0225 

Cimbala,  et  al.  (1981) 

u' 

156 

0.0167 

Present  Results 

u' 

500 

0.0284 

Present  Results 

v’ 

500 

0.0305 

The  general  trend  Is  faster  decay  rates  "with  increasing  Reynolds  number.  This 
is  especially  apparent  when  we  only  consider  data  from  one  facility  (e  g.,  the 
three  rates  of  Desruelle,  or  our  present  results).  There  is  much  scatter  in  the 
data  from  different  facilities;  the  decay  rate  is  probably  sensitive  to  background 
disturbances  in  the  freestreain. 

5.2.  Downstream  Structure 

As  downstream  distance  in  the  wake  is  increased,  fluctuations  appear  at  fre¬ 
quencies  lower  than  the  shedding  frequency.  Figure  5.3  shows  amplitude  spec¬ 
tra  at  several  downstream  positions  for  Re  =  150.  As  in  Figure  5.1,  the  spectra 
are  averaged  over  the  two  y -locations  where  the  slope  of  the  mean  velocity 
profile  is  a  maximum  (inflection  points).  We  observe  in  Figure  5,3  a  shift  to 
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lower  and  lower  frequencies  as  x/d  is  increased.  In  particular,  the  Karman  fre¬ 
quency  (/*  88  186  Hz)  Is  initially  the  only  significant  frequency.  At  xM  =  100  in 
Figure  5.3  a  broad  band  of  frequencies  has  arisen  with  the  most-prominent  peak 
at  /  «  90  Hz,  and  other  peaks  at  /  «  120,  83.  and  70  Hz.  By  x/d  =  200.  distur¬ 
bances  at  /  «  90  and  70  Hz  have  amplified  significantly,  as  well  as  the  subhar¬ 
monic  (/  =  83  Hz);  the  peaks  at  /  =  fk  =  166  Hz  and  at  /  =  120  Hz  have 
decayed.  With  further  increase  in  x/d,  the  band  of  prominent  disturbances 
shifts  to  lower  frequencies  and  broadens,  kt  x/d  -  400  the  most-prominent 
band  of  frequencies  is  centered  around  70  Hz,  while  all  frequencies  above  100  Hz 
have  disappeared.  By  x/d  =  750,  there  are  no  longer  any  outstanding  spikes; 
instead,  a  very  broad  hump  is  seen  for  6  <  /  <  75  Hz.  In  general,  following  the 
life-cycle  of  a  particular  frequency,  we  notice  that  the  fluctuation  amplitude  at 
that  frequency  first  rises,  then  decays;  the  lower  the  frequency,  the  further 
downstream  is  its  life-cycle.  In  Figure  5.4  we  have  plotted  spectral  amplitudes  at 
several  discrete  frequencies  versus  downstream  distance.  The  growth-decay 
cycle  is  clear. 

Notice  in  Figure  5.4  that  there  is  nothing  outstanding  about  the  subhar¬ 
monic,  /  =  83  Hz.  That  is  to  say.  disturbances  at  the  subharmonic  experience 
the  same  kind  of  growth  and  decay  as  at  any  other  nearby  frequency.  For  exam¬ 
ple,  the  growth-decay  cycle  at  /  «  90  Hz  leads  the  one  at  the  subharmonic  by 
about  20  diameters,  but  is  almost  identical  in  magnitude  and  shape. 

Figure  5.5  shows  amplitude  spectra  at  a  lower  Reynolds  number  of  140.  Here 
the  shift  to  lower  frequencies  is  more  "choppy;”  i.e.,  the  peaks  in  the  spectra  are 
not  as  broad  as  those  for  Re  =  150  -  it  appears  that  a  large  percentage  of  the 
fluctuation  energy  is  concentrated  at  discrete  frequencies  which  are  almost 
integral  multiples  of  25  Hz.  The  reason  for  this  behavior  is  not  entirely  clear.  As 
seen  in  Figure  B.8  (Appendix  B),  tunnel  disturbances  were  strong  at  25  Hz  and 


its  harmonics  relative  to  other  freestream  disturbances  -  the  wake  may  have 
"locked  in'*  at  these  frequencies  in  response  to  the  external  forcing.  Figure  5.6 
shows  spectral  amplitudes  at  several  of  these  discrete  frequencies  as  a  function 
of  xM  .  Their  growth-decay  cycles  are  very  similar  to  those  of  the  wake  at  Re  = 
150. 

Figure  5.7a)  also  shows  amplitude  spectra,  but  at  a  cylinder  Reynolds  number 
of  500,  where  the  wake  is  turbulent.  In  this  case  spectral  peaks  are  much  less 
prominent  because  of  fluctuations  at  all  frequencies.  A  broad  hump  is  discerni¬ 
ble,  however,  and  shifts  to  lower  frequencies  much  more  smoothly  than  for  the 
laminar  cases  above.  Notice  that  the  magnitude  of  the  entire  spectrum  decays 
with  downstream  distance,  as  energy  is  dissipated  by  turbulent  mixing.  Figure 
5.8a)  shows  growth-decay  cycles  at  Re  =  500;  we  don’t  really  see  regions  of 
growth,  but  rather  decay  at  each  discrete  frequency.  In  this  case  it  is  useful  to 
display  normalized  spectra  where  v'  is  divided  by  the  local  mean  velocity  defect 
A  U,  where  A 17  =  £/„  -  U0  =  UmW0,  instead  of  by  Um.  Figures  5.7b)  and  5.8b)  are 
normalized  versions  of  Figures  5.7a)  and  5.Ba)  respectively.  These  normalized 
plots  illustrate  the  growth-decay  cycles  much  more  clearly. 

In  Figure  5.7c)  we  show  normalized  spectra  at  all  twelve  downstream  posi¬ 
tions.  Here  we  have  normalized  frequency  as  well  as  spectral  amplitude,  in  the 
manner  discussed  in  Section  5.4  below.  The  fully  normalized  spectra  show  a 
self -similar  behavior  beyond  about  100  diameters.  The  center  of  the  spectral 
peak  in  these  coordinates  occurs  at  approximately  0  =  1.5.  The  band  of  prom¬ 
inent  frequencies  is  very  wide,  much  more  so  than  for  the  case  of  a  plane  mixing 
layer.  This  is  most  likely  due  to  the  much  slower  divergence  of  wakes  them  mix¬ 
ing  layers.  In  Section  5.4  we  compare  these  results  with  predictions  of  linear 
stability  theory. 


► 
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5.3.  The  Group  Phenomenon  ( 

By  studying  time-traces  of  velocity  fluctuations,  as  in  Figure  5.9,  at  Re  =  150, 
it  becomes  evident  that  the  secondary  structure  occurs  in  groups.  (Groups  are 
defined  as  packets  of  large  fluctuations,  lasting  for  several  cycles,  and  followed 

I 

by  periods  of  relative  calm.)  Recall  from  Chapter  1  that  Townsend  (1979)  had 

observed  this  phenomenon  for  the  turbulent  cylinder  wake.  His  results  show 

that  the  frequency  or  scale  of  the  fluctuations  is  constant  within  a  particular 

group  but  varies  from  group  to  group.  Here  we  wish  to  quantify  this  variation  -  * 

i.e.,  how  much  of  a  spread  in  frequency  is  there? 

Long-time-averaged  power  spectra  tell  us  little  about  the  instantaneous 
structure  in  a  flow  field.  However,  much  can  be  learned  by  examining  the  short- 
time-averaged  power  spectrum  of  each  group  of  structures.  Details  of  the 
numerical  technique  used  here  are  presented  in  Appendix  B.  Basically,  the  most- 
energetic  frequency  of  each  group  is  recorded;  a  probability  density  function  * 

(PDF)  of  the  spectral  peaks  is  then  calculated  for  all  these  frequencies,  weighted 
by  the  amplitude  of  the  fluctuations. 

Figure  5.10  shows  spectral-peak  probability  density  vs.  frequency  at  several  * 

downstream  positions  in  the  cylinder  wake  at  Re  =  150.  Close  to  the  cylinder,  the 
shedding  frequency  ( fk  «  166  Hz)  is  dominant  and  the  width  of  the  PDF  is 
extremely  small,  as  expected.  At  xM  =  100,  the  dominant  frequency  shifts  to  I 

about  90  Hz.  As  we  continue  downstream  we  see  a  shift  to  lower  frequencies  as 
well  as  further  broadening  of  the  frequency  band-width,  just  as  had  been 
observed  in  the  spectra  of  Figure  5.3.  ^ 

The  width  of  the  spectral-peak  PDF  is  significant  For  a  given  frequency  band 
to  have  non-zero  probability  density,  the  passage  frequency  of  at  least  one 
group  of  deflneable  structures  has  to  lie  within  that  frequency  band.  Whereas  9 

power  spectra  only  gi’  e  us  a  measure  of  the  energy  content  of  fluctuations  at  a 


9 
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given  frequency,  probability  densities  tell  us  that  identifiable  groups  of  energetic 
fluctuations  occur  at  discrete  frequencies.  The  broader  the  width  of  the  PDF.  the 
more  distributed  are  the  frequencies  of  these  groups.  Thus  our  results  support 
Townsend's  (1979)  observation  that  far-wake  structure  occurs  in  groups;  the 
frequency  within  a  group  being  fairly  constant,  but  varying  from  one  group  to 
the  next.  The  significance  of  this  result  will  be  more  apparent  in  the  following 
section,  where  we  compare  our  data  to  linear  stability  theory. 


Another  interesting  result  is  the  average  number  of  cycles,  rig,  within  a  group. 
The  original  Karman  street  of  course  is  practically  continuous,  and  rig  is 
extremely  large.  Figure  5.11  shows  how  n,  varies  with  downstream  distance 
beyond  100  diameters,  after  the  Karman  street  has  largely  decayed.  In  general. 
rig  lies  between  about  2  and  12,  but  seems  to  decrease  slowly  with  x/d . 


5.4.  Comparison  with  linear  Stability  Theory 

Inviscid  linear  stability  of  a  parallel  wake  profile  has  been  analyzed  by  many 
investigators,  both  in  the  temporal  and  spatial  modes.  Our  intent  here  is  to  com¬ 
pare  predictions  from  these  analyses  with  our  experiment  results.  Namely,  for 
the  measured  mean  velocity  profiles,  how  well  can  linear  theory  predict  which 
frequencies  are  most  prominent  in  the  wake?  It  is  convenient  to  express  the 
frequency  non-dimensionally  as  in  Figure  5.7c).  We  define 


Um  ' 


(5.1) 


where  /  is  the  dimensional  frequency,  and  6  is  the  wake  half-width,  defined  by 
the  y-location  where  U *  =  0.5.  U *  is  the  normalized  mean  velocity, 
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cr  = 

u0-  um 


(5.2) 


and  Uo  is  the  mean  centerline  velocity.  C/0  and  the  wake  hall-width  6  are  calcu¬ 
lated  by  fitting  a  Gaussian  profile  to  the  measured  mean  velocities,  as  described  * 

in  Appendix  B. 

For  comparison  with  experiment,  spatial  stability  analysis  is  preferred  over 
temporal  analysis,  as  has  been  pointed  out  by  Mattingly  and  Criminale  (1972). 

Kubota  (1983)  has  generated  numerical  solutions  for  both  cases;  inviscid  distur¬ 
bances  are  superposed  on  a  parallel  Gaussian  wake  profile.  Details  of  this  effort 
are  given  in  Appendix  C.  Figure  5.12  shows  spatial  growth  rate  as  a  function  * 

of  dimensionless  frequency  p  ,  and  for  several  values  of  W0.  the  normalized 
centerline  velocity  defect  in  the  wake.  (W0  =  1-U0/Um.)  Recall  that  for  spatial 
stability  analysis,  frequency  p  is  read,  while  wave  number  a  may  be  complex.  Fig-  • 

ure  5.13  shows  a  similar  family  of  curves  for  the  temporal  mode,  where  a  is  real 
and  (3  may  be  complex.  We  plot  spatial  growth  rate  g.  however,  using  the 
transformation  discussed  in  Appendix  C.  Notice  that  there  is  not  a  great  ^ 

discrepancy  between  the  two  modes.  We  are  not  concerned  here  with  predicting 
fluctuation  amplitudes  exactly;  our  main  concern  is  the  relative  amplitude  of 
fluctuations  at  various  frequencies.  Toward  this  end  either  mode  is  applicable; 
we  have  chosen  the  temporal  mode  for  convenience  in  our  analysis  below.  ® 

Another  important  thing  to  note  in  either  figure  is  that  as  velocity  defect 

decreases,  the  growth  curve  shifts  toward  higher  p.  Thus,  as  we  travel  down¬ 
stream  in  the  wake,  W0  decreases,  and  the  relevant  growth  curve  is  continually  • 

shifting  to  the  right,  as  well  as  decaying  in  magnitude. 


Kubota  (1983)  has  argued  that  the  growth  curves  shown  in  Figure  5.12  or  5.13 
may  be  extended  to  negative  growth  (de-amplification)  beyond  the  neutral  point, 


f3n.  This  is  extremely  important  in  the  following  analysis,  where  such  a  de- 
amplification  region  is  required  in  order  to  predict  growth-decay  cycles  of  large- 
scale  structures.  For  lack  of  details  in  this  region,  we  have  extrapolated  the 
growth  curves  in  Figure  5.13  by  straight  lines,  whose  slopes  are  matched  with 
those  at  the  neutral  points. 

Let  us  consider  which  prominent  frequency  fp  we  expect  to  find  in  our  meas¬ 
urements;  i.e.  we  wish  to  predict  fp  as  a  function  of  downstream  distance  * .  At 
first  glance  one  may  expect  Jp  to  equal  the  locally  most-amplified  frequency 
(the  peak  in  the  growth  curve).  This  turns  out  not  to  be  the  case  however,  due 
to  non-parallelism  of  the  mean  wake.  It  must  be  remembered  that  the  growth 
rates  shown  in  Figure  5.13  are  calculated  for  a  parallel  wake  (no  x -dependence); 
the  real  wake  of  course  is  not  parallel,  but  widens  downstream.  Obviously  the 
best  predictions  would  come  from  fully-non-parallel  stability  analysis,  where  one 
includes  downstream-widening  of  the  mean  wake;  nevertheless  we  expect  the 
curves  in  Figure  5.13  to  be  adequate  approximations  locally.  Hence,  a  scheme 
was  developed  to  approximate  non-parallelism  (  z -dependence)  using  only  paral¬ 
lel  stability  calculations.  Namely,  in  order  to  predict  how  large  a  disturbance  to 
expect  at  some  particular  frequency,  we  integrate  growth  rate  at  that  frequency 
with  respect  to  downstream  distance  z .  Doing  this  for  a  number  of  frequencies 
gives  us  disturbance  amplitude  as  a  function  of  frequency  for  each  downstream 
distance.  It  is  then  a  simple  matter  to  find  that  frequency  which  is  predicted  to 
be  the  most  prominent,  /,,,  as  a  function  of  z,  (i.e.,  the  frequency  whose  ampli¬ 
tude  is  greatest  at  a  given  z -location).  The  entire  process  is  outlined  below: 

1.  Measure  mean  velocity  profiles  at  several  downstream  positions  in  the  wake. 

2.  Calculate  (with  least-squares  fitting  to  a  Gaussian  wake  profile)  the  half¬ 
width  5  and  centerline  velocity  defect  W0  for  each  z -position.  (Examples  are 


shown  In  Figure  5.14.) 


3.  Fit  smooth  continuous  curves  through  W0  and  6  so  as  to  obtain  W0(x)  and 
6(x).  Examples  are  shown  in  Figures  5.15  and  5.16,  respectively. 

4.  The  appropriate  growth-rate  curve  g{x,f),  such  as  those  shown  in  Figure 
5.13,  is  chosen  locally  for  each  x  according  to  W0(x),  $(*).  and  Equation 
(5.1).  (  g  is  extrapolated  to  negative  growth  rates,  as  discussed  above.) 

5.  Knowing  g(x,f)  ,  we  assume  a  locally  exponential  growth  of  amplitude  A; 
namely, 

=  g(xJ)MxJ).  (5.3) 


We  can  predict  the  amplitude  of  fluctuations  at  frequency  /  and  at  location 
x  by  integrating  the  exponential  growth  from  x  =  0  to  z.  i.e., 


A(x.f)  =  A0  exp 


f  g(x\f)  dx • 


(5.4) 


where  A0(f )  is  the  initial  perturbation  amplitude  at  frequency  /,  and 
x*  =  x/5(x)  is  the  non-dimensionalized  downstream  distance.  In  Figure  5.17 
we  show  A(x,f)  versus  xM  for  a  number  of  fixed  frequencies  /  at  Re  =  150. 
We  have  assumed,  for  simplicity,  that  j40(/)  =  1  =  constant  for  all  /.  Notice 
the  growth-decay  cycles,  similar  to  our  experimental  results. 

8.  To  find  the  frequency  fp{x)  where  the  amplitude  given  by  Equation  (5.4)  is  a 
maximum,  we  simply  take  the  partial  derivative  of  A(x,f)  and  set  it  to  zero: 


or.  at  /  =  /,  . 

/  d9(<Hdx*  =  0.  (5.5) 

7.  Compare  this  predicted  peak  frequency,  fp(x),  to  the  measured  peak  fre¬ 
quency  at  each  x-location. 

Let  us  now  examine  our  experimental  results  at  Re  =  150  in  light  of  the  above 
discussion.  We  show  in  Figure  5.18  a  plot  of  frequency  versus  xM.  The  circles 
represent  experimental  frequencies  at  which  the  most-prominent  disturbances 
were  measured.  The  vertical  lines  represent  70%  confidence  levels,  and  give  the 
reader  a  feel  for  the  band-width  of  prominent  frequencies.  These  observed  fp 
should  be  compared  to  the  predicted  fp  as  described  above  —  given  here  as  the 
solid  line.  Also  shown  in  Figure  5.18  are  the  locally  neutral  frequency  (dashed 
line),  and  the  frequency  corresponding  to  the  locally  most-amplified  disturbance 
(dotted  line).  Note  that,  because  of  the  integration,  the  frequency  which  is 
locally  most-amplified  is  not  necessarily  the  frequency  with  the  largest 
(integrated)  amplitude. 

Close  to  the  cylinder,  the  Karman  frequency  is  the  most-prominent  observed 
frequency,  as  expected.  Comparison  with  linearized  theory  is  futile  in  the  very 
near  wake  where  the  amplitude  of  fluctuations  can  exceed  30%  (Nishioka  and 
Sato  (1978)).  Note  that  the  Karman  vortices  are  initially  formed  by  a  mechan¬ 
ism  much  different  than  far-wake  instability.  The  mechanism  consists  of  a 
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complex  combination  of  unsteady  boundary-layer  separation  and  near-wake  ins¬ 
tability,  and  is  not  yet  completely  understood.  We  therefore  do  not  expect  fk  to 
match  our  far-wake  predictions.  Indeed,  it  is  interesting  in  Figure  5. IB  that  the 
Karman  frequency  is  well  above  our  predicted  peak  frequency.  In  fact,  fk  even 
lies  above  the  locally  neutral  frequency  and  is  thus  in  the  region  of  damping, 
according  to  stability  theory.  This  is  consistent  with  our  observation  that 
fluctuations  at  Karman  frequency  /*  decay  exponentially  in  the  region  25  <  x/d 
<  125. 

As  disturbances  at  fk  decay  rapidly,  lower  frequencies  take  over;  beyond  x/d 
-  100  the  agreement  between  predicted  and  observed  frequencies  is  quite  good. 

Figure  5.19  shows  a  similar  comparison  for  Re  =  140.  Again,  fk  (»  149  Hz) 
lies  above  the  neutral  curve;  fluctuations  at  fk  decay  beyond  our  first  measure¬ 
ment  station  ( x/d  =  40).  Lower-frequency  disturbances  become  prominent 
beyond  100  diameters,  and  the  agreement  with  the  predicted  frequencies  is 
again  quite  good.  Subharmonic  fk/Z  «  75  Hz  seems  to  lock  in  for 
125  <  x/d  <  250.  Beyond  x/d  =  250,  fk/Z  lies  above  the  neutral  curve,  lower- 
frequency  fluctuations  arise  as  the  subharmonic  decays.  The  70%  confidence 
level  bars  are  quite  short;  this  implies  very  narrow  band-widths,  as  had  been 
observed  previously  in  the  spectra  at  this  Reynolds  number.  (See  Figure  5.5.) 

Figure  5.20  shows  the  comparison  between  theory  and  experiment  for  the 
case  of  a  turbulent  cylinder  wake  at  Re  =  500.  Here,  fk  is  approximately  630  Hz, 
and  is  ofl  the  scale.  In  this  case  fk  is  prominent  only  in  the  very  near  wake. 
Already  by  xM  =  50.  a  lower-frequency  band  (centered  around  fp  »  290  Hz)  is 
prominent.  Notice  that  the  observed  fp  decreases  quite  smoothly  with  x  as  we 
travel  downstream.  90%  confidence  levels  have  been  chosen;  the  band  of  prom¬ 
inent  frequencies  is  very  broad.  Interestingly,  our  observed  Jp  lies  above  the 
predicted  fp  ;  in  fact,  it  follows  almost  exactly  the  locally  neutral  frequency  /„. 


It  is  unclear  why  our  prediction  scheme  works  so  well  for  Re  =  140  and  150,  but 
not  for  Re  =  500.  At  Re  =  500  the  wake  is  turbulent,  and  fluctuations  are  strong 
at  all  frequencies  in  the  range  of  interest.  Recall  from  Figures  5.7a)  and  5.8a) 
that  fluctuations  at  any  given  frequency  never  really  have  a  chance  to  "grow," 
since  the  entire  power  spectrum  decays  with  xM-  The  most-prominent  fre¬ 
quency  fp  is  then  the  one  which  suffers  the  least  integrated  decay  as  the  distur¬ 
bances  travel  downstream.  That  frequency  is  apparently  the  locally  neutral  fre¬ 
quency /„.  as  Figure  5.20  shows.  No  attempt  has  been  made  in  our  calculations 
to  account  for  the  additional  rapid  decay  of  fluctuations  due  to  turbulent 
stresses.  Perhaps  a  more  rigorous  model  is  necessary.  We  may  point  out  here 
that  similar  results  have  been  shown  recently  by  Wygnanski,  et  al.  (1983)  for  the 
case  of  a  plane  turbulent  wake  at  Re  =  4000.  Namely,  the  observed  prominent 
frequency  matches  well  with  the  neutrally  amplified  frequency,  as  predicted  by 
linear  stability  theory.  Ho  and  Huerre  (1984)  have  observed  this  agreement  for 
the  case  of  a  plane  turbulent  mixing  layer  as  well.  Lessen  and  Singh  (1974), 
using  eddy  viscosity  models,  have  postulated  a  similar  behavior  for  axisym- 
metric  turbulent  jets  and  wakes. 

Indeed,  it  is  surprising  that  such  a  simple  prediction  scheme  as  employed 
here  (with  local  parallelism)  works  as  well  as  it  does.  We  may  conclude  from  this 
exercise  that  two-dimensional  inviscid  parallel  hydrodynamic  stability  theory, 
when  interpreted  properly,  is  adequate  to  predict  which  frequencies  are 
expected  in  a  far  wake.  Furthermore,  this  agreement  supports  our  hypothesis 
that  structure  appearing  far  downstream  of  a  bluff  body  is  the  result  of  hydro- 
dynamic  instability  of  the  mean  velocity  profile,  and  is  not  directly  dependent  on 
the  shed  vortices  of  the  Karman  vortex  street. 
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S.5.  A  Look  at  y-Dependence 

As  discussed  in  Chapter  2,  “on-line"  spectral  analysis  was  also  done  with  a 
single  inclined  hot  wire  and  a  digital  spectrum  analyzer.  In  particular,  we  are 
interested  in  the  dependence  of  amplitude  spectrum  on  distance  from  the  wake 
centerline  (y -direction).  Figure  5.21  shows  long-time-averaged  amplitude  spec¬ 
tra  in  a  circular  cylinder  wake  at  Re  =  300,  at  x/d  -  750,  and  at  several  y- 
locations.  Output  from  the  hot-wire  anemometer  was  input  directly  to  the  spec¬ 
trum  analyzer,  without  any  linearization;  the  amplitude  scale  is  therefore  in 
arbitrary  (log)  units.  (Velocity  fluctuations  are  very  small  this  far  downstream; 
hot-wire  response  is  assumed  to  be  linear.) 

Close  to  the  centerline  (0  <  y  M  <  15).  no  peak  in  the  spectrum  can  be 
found.  For  y/d>  15  however,  a  broad-band  peak  in  the  spectrum  exists;  furth¬ 
ermore,  the  band  shifts  to  lower  frequencies  and  becomes  narrower  as  we  travel 
outward  in  the  wake.  The  most-prominent  frequency  fp  is  plotted  as  a  function 
of  y  M  in  Figure  5.22;  an  almost  linear  decrease  in  fp  is  seen. 

We  can  explain  this  frequency  shift  as  follows:  The  wake  profile  at  a  particular 
z-location  is  unstable  to  a  broad  band  of  frequencies.  Fluctuations  at  high  fre¬ 
quencies  may  be  measured  occasionally,  the  scale  of  which  is  smaller  than  the 
local  wake  scale.  We  can  think  of  these  structures  as  originating  upstream  of 
our  measurement  position.  Because  of  their  smedler  scale,  they  will  be  sensed 
more  readily  by  a  probe  nearer  the  wake  centerline.  When  the  probe  is  near  the 
outer  edge  of  the  wake,  it  can  only  sense  structures  of  larger  scale  (lower  fre¬ 
quency),  which  can  be  thought  of  as  originating  further  downstream.  These 
ideas  are  sketched  in  Figure  5.23.  Sketch  a)  shows  high-frequency  structures 
which  originate  upstream  and  are  sensed  at  probe  position  y  ^  but  not  at  y2  , 
sketch  b)  shows  lower-frequency  structures  which  can  be  sensed  at  probe  posi¬ 
tion  yz 
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5.6.  Simultaneous  Smoke-Wire  and  Cross-Wire  Records 

The  microprocessor  was  programmed  such  that  a  smoke-wire  photograph 
could  be  taken  at  any  specified  time  during  the  data-collection  interval.  As  illus¬ 
trated  in  Figure  2.2,  the  microprocessor  first  sends  a  pulse  to  the  relay  unit, 
which  turns  on  the  smoke  and  opens  the  camera  shutter  (the  room  is  dark,  of 
course).  After  a  specified  delay  to  allow  the  smoke  to  convect  completely  into 
the  field  of  view,  the  strobe  is  triggered  to  flash  once.  Subsequently  the  smoke  is 
turned  off,  and  the  camera  shutter  closed. 

The  sequence  outlined  above  enables  us  to  display  smoke-wire  and  hot-wire 
records  simultaneously;  an  example  is  shown  in  Figure  5.24  for  a  cylinder  wake 
at  Re  =  150.  The  smoke-wire  is  at  x/d  =  275,  while  the  cross-wire  is  positioned  at 
xM  =  400  and  at  y  At  =  -4.  Note  that  in  order  to  avoid  damage  to  the  hot  wires, 
the  cross-wire  assembly  was  moved  slightly  out  of  the  plane  of  smoke.  The  time 
trace  in  Figure  5.24  has  been  displayed  ••backwards’’  (i.e.,  time  increasing  to  the 
left),  so  that  direct  comparison  between  the  flow-visualization  photograph  and 
time-trace  is  possible.  The  scales  have  been  matched  according  to 


where  Uc  is  the  local  convection  velocity  (assumed  to  be  the  mean  velocity  at  the 
location  of  the  probe).  In  this  example,  d  =  0.0016  m,  Um  =  1.52B  m/s.  and 
Uc  =  1.397  m/s.  The  time  axis  has  also  been  shifted  such  that  the  strobe  flash 
occurs  at  t  =  0. 

Because  of  the  integration  effect  of  streaklines,  as  discussed  in  Chapter  3,  we 
do  not  expect  the  time-trace  to  exactly  match  the  distortions  of  streaklines  in 
the  photograph.  Nevertheless,  a  general  agreement  is  certainly  present.  Notice 
for  example  the  burst  of  orderly  structure  at  -  60  <  f  <  0;  large-scale  vortices 
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are  clearly  seen  in  the  photograph  at  400  <  x/d  <  460.  Just  ahead  of  the  probe 
(350  <  x/d  <  390)  there  appears  in  the  photograph  a  smaller-scale  group  of 
vortices.  This  is  also  recognizable  in  the  time-trace  as  a  very  weak  higher- 
frequency  oscillation.  Short-time  spectra  of  the  fluctuations  shown  in  Figure 
5.24  give  frequencies  of  64  Hz  for  the  larger-scale  structure,  and  84  Hz  for  the 
smaller-scale  structure,  Both  of  these  frequencies  are  within  the  range  of 
amplified  frequencies,  according  to  linear  stability  theory,  as  can  be  seen  in  Fig¬ 
ure  5.18  at  x/d  -  400. 

Both  the  time-trace  and  smoke-wire  photograph  of  Figure  5.24  support 
Townsend’s  suggestion  that  far-wake  structure  appears  as  groups  of  several  vor¬ 
tices.  Such  a  hypothesis  is  consistent  with  the  explanation  that  downstream 
structure  arises  because  of  wake  instability  of  the  mean  velocity  profile  and  not 
directly  because  of  vortex  amalgamation.  Furthermore,  Figure  5.24  shows  that 
far-wake  structure  convects  at  the  local  mean  velocity;  it  does  not  appear  to 
have  unequal  group  and  phase  velocities. 

5.7.  Porous  Flat  Plate  Wake  Surveys 

Cross-wire  surveys  were  also  taken  in  the  wake  of  a  47%  solid  flat  plate  at 
Re*  =  5000,  where  d  is  the  width  of  the  plate.  Mean  velocity  profiles  are  shown  in 
Figure  5.25.  As  expected,  very  near  to  the  plate  the  profile  is  a  sharp  “top  hat” 
which  gradually  relaxes  as  the  shear  layers  on  either  side  of  the  wake  begin  to 
merge.  By  ten  diameters  the  profiles  are  far-wake-like,  and  have  been  fitted  here 
with  Gaussians. 

Amplitude  spectra  of  v'/U„  are  plotted  in  Figure  5,26  for  several  z -locations; 
each  spectrum  is  averaged  at  the  two  inflection  points  of  the  mean  velocity 
profile,  as  discussed  previously.  At  x/d  =  1,  spectral  energy  at  low  frequencies  is 
small,  but  there  is  a  noticeable  bump  centered  around  410  Hz.  Recall  from  our 
flow-visualization  results  that  the  plate  does  not  appear  to  shed  Karman-type 
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vortices;  the  peak  in  the  spectrum  may  instead  be  associated  with  the  plane 
shear  layers  which  exist  initially  on  the  outer  edges  of  the  wake.  We  can  calcu¬ 
late  the  local  Strouhal  number  of  one  of  these  shear  layers: 

St  =  (5.7) 


where  9sl  is  the  shear  layer  integral  thickness, 


6sl  -  J 


'  u  -  ux' 

to 

1 

to 

U2  -  C/i 

1  Ua-u ,  J 

dy. 


(5.8) 


and  V  is  the  average  of  Ux  and  Ue.  At  x/d  =  1  in  Figure  5.26  we  have  Ux  -  2.39 
m/s.  Uz  =  4.54  m/s.  and  eSi  =  0.0507  cm.  For  /  =  410  Hz,  Equation  (5.7)  gives 
St  =  0.06.  which  is  consistent  with  the  observations  of  Ho  and  Huerre  (1984)  for 
a  growing  free  shear  layer.  Namely,  they  report  that  St  =  0.032  when  the  shear 
layer  first  begins  to  oscillate,  and  St  =  0.079  further  downstream  as  the  shear 
layer  grows.  Our  value  lies  between  these  two  extremes,  which  confirms  our 
speculation  that  the  peak  at  /  =  410  Hz  is  associated  with  the  shear  layers  on 
either  side  of  the  wake,  and  not  with  any  kind  of  Karman-type  vortex  shedding. 

At  four  diameters  downstream  the  amplitude  spectrum  is  devoid  of  any  out¬ 
standing  peaks,  except  for  a  very  broad  band  centered  at  approximately  150  Hz; 
the  peak  at  410  Hz  has  disappeared.  We  shall  designate  the  range  4  <  x/d  <  10 
as  a  transition  range,  where  the  stability  problem  changes  from  that  of  two 
(independent?)  plane  shear  layers  to  that  of  a  plane  far  wake. 

At  x/d  =10  the  spectrum  has  increased  in  magnitude  at  the  lower  frequen¬ 
cies,  and  in  particular  there  is  a  broad  hump  centered  around  75  Hz.  Beyond  10 
diameters,  in  the  “far  wake"  region,  the  peak  in  the  spectrum  shifts  to  lower 
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frequencies  (larger  scales)  as  downstream  distance  increases. 

Figure  5.2?  shows  a  comparison  between  observed  and  predicted  prominent 
frequencies  fp  ;  the  predicted  values  come  from  our  locally-parallel  stability 
analysis,  as  discussed  in  Section  5.4.  The  agreement  is  quite  good  beyond  xAL~ 
10,  in  th§  far-wake  region.  75%  confidence  levels  are  also  plotted  in  Figure  5.2? 
to  indicate  the  band-width  of  prominent  frequencies. 

Direct  comparison  between  circular  cylinder  wakes  and  wakes  of  other  bluff 
bodies  can  be  accomplished  by  equating  Reynolds  number  based  on  momentum 
thickness  e,  where 


*=/ 


(5.9) 


For  the  case  of  the  47%  solid  plate  discussed  here,  we  have  d  =  1.7B  cm,  6Ai  =  • 

0.35,  and  Um  -  4.54  m/s;  Re*  is  therefore  equal  to  1750.  For  our  circular 
cylinder  wakes  we  found  e/d  »  0.8  ;  Re*  is  thus  60%  of  Re* . 
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Chapter  6 

Three-Dimensional  Effects 

I 

6.1.  Unforced  Three-Dimensionality 

Up  to  this  point  our  discussion  has  been  confined  to  the  x-y  plane;  now  let  us 
consider  three-dimensionality.  Possibly  the  earliest  experiments  on  far-wake  I 

structure  are  those  of  Townsend  (1956)  and  of  Grant  (1958).  Grant’s  work,  an 
extension  of  Townsend's,  consists  of  measured  velocity  correlations  in  the  far 
wake  («  500  diameters  )  of  a  two-dimensional  cylinder  at  Re  =  1300.  A  pecu¬ 
liar  thraa-tHmansional  structure  was  inferred  from  these  time-averaged  meas¬ 
urements,  and  was  christened  the  "vortex-pair  eddy,”  later  called  the  "double- 

I 

roller  eddy"  by  Townsend  (1970).  More  detailed  calculations  by  Payne  and  Lum- 

ley  (1967)  of  Grant’s  data  yielded  a  similar  structure.  * 

Since  that  time,  neither  the  instantaneous  nature  of  this  structure  nor  its 
origin  have  been  adequately  explained,  although  there  has  been  some  specula- 

I 

tion  (see  for  example.  Keffer  (1965),  Townsend  (1970),  Roshko  (1976)  and  Town¬ 
send  (1979)).  In  particular,  Roshko  (1978)  suggested  the  structure  may  actually 
be  the  time-averaged  superposition  of  vortex  loops,  formed  by  the  pinching  off 
and  joining  together  of  vortices  from  opposite  sides  of  the  street.  That  long-  • 

time-averaging  seriously  distorts  one’s  interpretation  of  the  instantaneous 
structure  was  also  mentioned  by  Townsend  (1979).  He  suggested  that  time¬ 
averaging  of  the  large-scale  velocity  patterns  “makes  them  appear  more  com-  I 

plex  because  of  the  superposition  of  patterns  from  eddies  at  all  stages  of  ...  [his 
proposed]  growth-decay-renewal  cycle.”  An  excellent  discussion  of  some  of 
these  ideas  is  offered  byYflezien  (1981). 


Most  recently  Mumford  (1983),  using  a  pattern-recognition  technique,  sug¬ 
gests  that  the  double-roller  eddies  may  be  even  more  complex  than  previously 
thought.  He  found  that  the  structures  are  often  confined  to  one  side  of  the 
wake  centerplane,  rather  than  extending  across  the  entire  wake.  Furthermore, 
eddies  of  similar  type  tend  to  occur  in  groups  of  two  or  more,  one  after  another 
in  the  streamwise  direction 

Regardless  of  its  exact  nature,  it  is  interesting  that  the  characteristic  large 
structure  far  downstream  in  a  plane  wake  appears  to  be  three-dimensional. 

The  smoke-wire  flow-visualization  setup  was  modified  for  plan  views  by  align¬ 
ing  the  cylinder  parallel  to  the  smoke-wire.  Figure  6.1  shows  both  edge  (x-y 
plane)  and  plan  (x-x  plane)  views  of  the  cylinder  wake  at  Re  =  140.  (The  two 
views  were  recorded  for  the  same  tunnel  conditions,  though  not  simultane¬ 
ously.)  The  edge  view  illustrates  the  decay  of  the  Karman  vortex  street  and  the 
subsequent  growth  of  secondary  structure.  The  smoke-wire  was  positioned  at 
xM  =  8.  100,  and  200  in  the  manner  discussed  previously. 

For  the  plan  view,  the  smoke-wire  was  at  the  same  three  x -positions  as  above, 
but  at  yM  =0  (i.e.,  in  the  center-plane  of  the  wake),  and  parallel  to  the 
cylinder.  In  photo  a)  one  immediately  notices  the  skewed  angle  at  which  vortices 
are  shed  from  the  body.  The  reason  for  this  is  not  entirely  clear,  but  similar 
observations  have  been  reported  by  many  others.  (See  for  example  Tritton 
(1969),  Gerrard  (1966),  Nishioka  and  Sato  (1978).)  Our  own  experiments  indi¬ 
cate  the  slantwise  shedding  to  be  quasi-stable;  i.e.,  vortices  are  sometimes  shed 
parallel,  sometimes  slanted  one  way,  and  other  times  the  opposite  way.  Subtle 
non-uniformities  of  the  freestream,  end  conditions,  or  the  body  itself  are  the 
most  likely  candidates  for  triggering  transitions  between  the  three  "modes.”  Fig¬ 
ure  6.2  shows  an  unusual  case  where  we  have  vortices  shed  slantwise  at  two 
opposing  angles  with  an  "elbow"  in  between. 
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At  any  rate,  in  Figure  6.1  the  shed  vortices  are  parallel  to  each  other,  with  no 
appreciable  spanwise  structure.  After  the  decay  ot  the  primary  street  however, 
three-dimensionality  is  seen  in  the  secondary  structure  beyond  100  diameters 
(photos  b  and  c).  The  honeycomb-like  pattern  also  is  skewed  in  these  photo¬ 
graphs;  but,  just  as  with  the  shed  vortices,  the  skewness  comes  and  goes.  The 
size  of  the  cells  is  approximately  20  cylinder  diameters. 

What  exactly  is  this  structure?  What  is  the  mechanism  by  which  it  develops? 
How,  if  at  all,  does  it  relate  to  Grant’s  “vortex-pair  eddy?" 

Figure  6.3  may  provide  a  hint  for  at  least  the  second  of  these  questions.  For 
Re  =  150,  with  the  smoke-wire  at  xM  -  100,  we  see  the  gradual  formation  of 
three-dimensional  structure.  By  moving  the  smoke-wire  out  of  the  wake  center- 
plane  (  yAt  =  2),  we  are  able  to  visualize  vortex  lines  from  just  one  side  of  the 
secondary  street.  In  this  particular  case,  the  vortices  are  initially  straight  and 
parallel.  A  waviness  quickly  develops,  with  successive  vortices  16 0°  out  of  phase. 
The  amplitude  of  this  perturbation  grows  downstream,  eventually  to  where  one 
vortex  overlaps  the  next.  Presumably  the  other  side  of  the  street  has  a  similar 
experience;  the  interaction  between  the  two  sides  may  be  quite  complex,  but  it  is 
not  inconceivable  that  such  an  interaction  could  culminate  in  vortex  "loops.” 
As  was  suggested  by  Roshko  (1976),  iong-time-averaged  correlations  of  these 
loops  could  be  interpreted  as  Grant’s  "vortex-pair  eddies.”  Some  evidence  that 
vortex  loops  con  exist  in  a  wake  is  provided  by  Breidenthal  (1980).  A  three- 
dimensionally-perturbed  splitter  plate  in  a  shear-layer  facility  was  used  to  pro¬ 
duce  the  loops,  which  persist  for  large  downstream  distances. 

For  higher  Reynolds  numbers,  where  the  wake  becomes  turbulent,  smoke 
visualization  is  encumbered  by  small-scale  structure.  In  Figure  6.4  however,  at 
Re  =  190,  we  are  still  able  to  recognize  wavy  structures  similar  to  those  of  the 
laminar  wake. 
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8.2.  Forcing  Attempts 

We  should  point  out  here  that  Figures  6.1  and  6.3  represent  exceptional  '' 

cases.  Only  rarely  is  the  three-dimensionality  as  strong  as  in  these  photographs. 
Nevertheless,  it  seems  that  there  exists  an  Inherent  three-dimensional  instabil¬ 
ity,  which  under  proper  conditions  can  become  quite  significant.  We  made  a  » 

number  of  attempts  at  enhancing  three-dimensionality;  for  the  most  part  these 
attempts  were  not  successful. 

The  first  such  attempt  was  to  modify  the  cylinder  itself.  Paint  was  sprayed  on  * 

the  cylinder  at  regular  intervals,  causing  periodic  variation  of  its  diameter.  The 
hope  was  to  generate  fixed-wavelength  disturbances,  which  would  initiate  ord¬ 
erly  three-dimensional  structure  downstream.  Unfortunately  even  a  very  small  I 

diameter  perturbation  («  2%)  significantly  disturbs  the  vortex  shedding 
mechanism,  and  leads  to  very  irregular  unsteady  flow  patterns  downstream. 

Orderly  three-dimensionality  Is  not  enhanced. 

I 

Acoustic  excitation  was  also  attempted.  A  loudspeaker  was  flush-mounted  to 
the  tunnel  wall  at  about  150  diameters,  such  that  the  speaker  axis  was  parallel 
to  that  of  the  cylinder.  Theoretically,  for  certain  frequencies  of  excitation,  a  , 

standing  wave  pattern  can  be  generated  between  the  walls  of  the  tunnel,  thereby 
imposing  a  steady  fixed-wavelength  disturbance  on  the  wake.  Practically,  how¬ 
ever,  it  is  extremely  difficult  to  generate  a  steady  standing  wave;  therefore  this 
technique  was  also  unsuccessful  at  enhancing  three-dimensionality.  ' 
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Chapter  7 

Discussion  and  Conclusions 

7.1.  Interpretation  of  Streakline  Patterns 

We  found  smoke-wire  flow  visualization  to  be  very  effective,  as  well  as  easy-to- 
use.  in  the  study  of  wakes  of  bluff  bodies.  We  must  emphasize,  however,  that 
extreme  caution  should  be  taken  when  interpreting  streakline  patterns.  In  par¬ 
ticular,  when  the  classical  case  of  a  Karman  vortex  street  is  visualized  by  place¬ 
ment  of  smoke-wires  at  several  downstream  locations,  the  confusing  integration 
effect  of  streaklines  becomes  immediately  apparent.  The  conclusion  to  be 
drawn  from  this  exercise  is  that  in  an  unsteady  flow  it  is  best  to  introduce  flow 
tracers  at  not  just  one,  but  several  positions  in  the  flow.  Only  by  studying  a  col¬ 
lage  of  photographs  arranged  in  the  manner  described  in  this  paper  can  one  be 
confident  that  he  is  not  being  misled  by  the  integration  effect  of  streaklines. 

7.2.  Decay  of  the  Karman  Vortex  Street 

It  has  been  shown,  with  flow  visualization  as  well  as  with  hot-wire  measure¬ 
ments.  that  the  Karman  vortex  street  shed  in  the  wake  of  a  circular  cylinder 
does  not  persist  indefinitely,  but  rather  decays  exponentially  with  downstream 
distance.  In  the  present  experiments,  the  decay  is  so  rapid  that  fluctuating 
velocity  measurements  at  the  Karman  frequency  are  lost  in  background  "noise” 
by  100  to  150  diameters  at  Re  =  150.  The  decay  is  even  faster  for  turbulent 
wakes  where  small-scale  structure  is  prevalent. 

Such  rapid  decay  is  by  no  means  a  new  revelation;  however  its  full  impact  has 
not  been  realized  up  to  now  because  of  the  misinterpretation  of  streakline  flow 
visualization,  as  discussed  in  Chapter  3.  Reasons  for  the  decay  are  certainly  a) 


viscosity,  and  b)  de-amplification  in  the  sense  of  hydrodynamic  stability.  In  our 
numerical  streakline  simulations  only  the  first  of  these,  i.e.  viscosity,  is 
included.  The  decay  rate  inferred  from  Figure  3.2  is  not  as  rapid  as  that  for  the 
experimental  case  of  Figure  3.1;  this  shows  that  viscosity  is  not  the  only 
mechanism  of  decay.  Desruelle  (1983)  shows  evidence  that  the  second  mechan¬ 
ism,  hydrodynamic  de-amplification,  accounts  for  at  most  77%  of  the  vortex 
street  decay  rate.  Other  possibilites  contributing  to  the  decay  include  c)  cancel¬ 
lation  of  vorticity  from  opposite  sides  of  the  street,  and  d)  the  Karman  vortices 
may  be  located  differently  with  respect  to  the  centerline  than  a  regular  eigen¬ 
function.  This  final  point  was  suggested  by  Desruelle,  where  he  shows  that  as  6 
increases  with  x,  the  lateral  spacing  b  between  Karman  vortices  does  not 
"keep  up;"  Karman  vortices  that  initially  sit  on  the  edge  of  the  wake  move  closer 
to  its  center  as  the  wake  grows.  More  detailed  studies  of  this  phenomenon  and 
its  effect  on  the  decay  of  a  vortex  street  have  not  been  attempted  to  date. 

7.3.  Downstream  Structure 

As  the  wake  widens  with  downstream  distance,  the  scale  of  large  structures 
must  also  increase.  The  mechanism  by  which  this  is  accomplished  is  hydro- 
dynamic  instability  of  the  developing  mean  wake  profile.  Frequencies  are  selec¬ 
tively  amplified  and  then  damped,  according  to  the  local  growth  rates,  as  deter¬ 
mined  by  the  local  mean  wake  profile.  Observed  prominent  frequencies  agree 
fairly  well  with  predictions  of  parallel  inviscid  linear  stability  theory,  provided 
that  streamwise  growth  of  the  wake  is  taken  into  account. 

The  distribution  of  frequencies  which  are  prominent  broadens  with  down¬ 
stream  distance.  The  structure  appears  as  groups  of  several  vortices,  separated 
by  regions  of  random  fluctuations;  the  frequency  (or  scale)  of  vortices  within  a 
group  is  constant,  while  it  can  vary  considerably  from  group  to  group.  Such  a 
pattern  agrees  intuitively  with  what  one  expects  in  a  hydrodynamically  unstable 
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system.  As  a  disturbance  grows  It  convects  downstream  until  it  Is  no  longer 

> 

amplified,  but  damped.  Meanwhile,  in  a  somewhat  random  fashion,  disturbances 
at  other  frequencies  are  amplified,  perhaps  originating  at  different  x  -positions, 
and  go  through  similar  life-cycles. 

> 

The  above  scenario  is  appropriate  for  turbulent  as  well  as  laminar  wakes.  In 
the  turbulent  case,  additional  complications  vise  because  of  small-scale  tur¬ 
bulent  structure.  Townsend  (1966)  proposes  his  equilibrium  hypothesis  as  a 
mechanism  for  the  growth-decay  cycles  of  large  structures.  Namely,  large  * 

eddies  develop  during  a  period  of  quiescence,  when  the  turbulence  intensity  is 
low.  The  large  eddies  grow  in  strength  and  cause  rapid  entrainment,  which  leads 
to  an  Increase  of  turbulent  intensity  at  all  scales.  Turbulent  motion  of  smaller  l 

scales  resists  the  growth  of  large  eddies  by  absorbing  some  of  their  energy.  The 
large  eddies  therefore  decay,  and  another  period  of  quiescence  begins.  Further 
downstream,  where  the  wake  is  wider,  larger-scale  (lower-frequency)  structure  ' 

emerges,  and  the  cycle  recurs.  Our  results  are  not  inconsistent  with  Townsend’s 
equilibrium  hypothesis.  We  stress  in  addition  the  important  role  of  local  hydro- 
dynamic  stability  in  the  growth  and  decay  of  large  structures. 

> 

Let  us  now  address  the  phenomenon  of  vortex  amalgamation.  As  discussed  in 
the  introduction,  there  seems  to  be  some  dispute  as  to  the  role  of  vortex  pairing 
in  the  changeover  to  larger-scale  structure  in  wakes.  There  are  those  who  argue 
that  pairing  instability  is  paramount  in  this  changeover.  It  is  our  contention 
however  that  vortex  amalgamation  is  not  necessary  for  the  growth  of  larger 
structure;  but  rather,  hydrodynamic  instability  of  the  developing  mean  velocity 
profile  results  in  the  selective  amplification  of  frequencies  which  are  locally  * 

unstable.  If  amalgamation  does  occur  it  Is  purely  incidental,  not  the  driving 
mechanism  of  the  flow.  Wakes  of  porous  flat  plates  are  particularly  supportive 
of  this  statement.  Wake  instability  alone  establishes  a  street-like  structure  » 


downstream  -  no  shed  vortices  even  exist. 

Nevertheless,  we  do  observe  the  subharmonic  fk/2  to  be  a  popular  frequency 
in  the  data,  particularly  at  125  <  xM  <  250  for  the  cylinder  wake  at  Re  =  140. 
At  x/d  «  100,  the  primary  (Karman)  street  has  not  completely  decayed;  it  is 
likely  that  disturbances  at  the  subharmonic  fk; 2  are  present  because  of  the 
residual  Karman  vortices.  If  such  is  the  case,  these  disturbances  have  a  “head- 
start”  at  being  amplified  when  compared  to  those  at  frequencies  not  close  to 
fk/2.  Thus  fk/2  disturbances  may  be  chosen  to  be  amplified,  even  though  they 
are  not  necessarily  the  maximally-amplified  disturbances  predicted  by  stability 
theory.  The  residual  Karman  vortices  will  pair,  not  due  to  a  pairing  instability  in 
and  of  itself,  but  merely  in  response  to  the  growing  subharmonic  disturbance. 
For  the  most  part,  the  frequency  of  downstream  fluctuations  shifts  gradually, 
particularly  for  turbulent  wakes.  At  Re  =  500,  for  example,  no  identifiable 
trends  exist  that  would  indicate  amalgamation. 

Indeed,  except  perhaps  for  the  special  case  of  Re  =  140,  vortex  pairing  was 
not  observed  by  this  experimenter.  At  Re  =  150  or  500,  for  example,  there  is 
nothing  outstanding  about  the  subharmonic  frequency.  That  is  to  say,  distur¬ 
bances  at  the  subharmonic  experience  the  same  kind  of  growth  and  decay  as  at 
any  other  nearby  frequency. 

A  plot  of  a2/al  (ratio  of  secondary-street  scale  to  that  of  the  primary 
street)  versus  Reynolds  number,  as  shown  for  example  by  Matsui  and  Okude 
(1980),  has  a  great  amount  of  scatter,  particularly  among  different  flow  facili¬ 
ties.  Furthermore,  a2/ai  is  exactly  equal  to  two  only  for  isolated  cases.  Hydro- 
dynamic  stability  theory  adequately  explains  the  scatter.  As  indicated  by  Figure 
5.13,  the  wake  at  a  given  location  is  unstable  to  a  broad  band  of  fequencies.  If 
there  exist  relatively  large  background  disturbances  at  specific  frequencies 
within  this  band,  these  will  most  likely  be  the  ones  selected  to  amplify.  It  is  then 


clear  why  such  facility-dependency  exists.  Moreover,  because  a  broad  band  of 
scales  are  amplified  in  the  wake,  and  because  the  band  of  amplified  frequencies 
shifts  to  the  left  withx,  az  can  not  be  defined  unambiguously;  a  plot  of  a2/o,  is 
thus  meaningless,  perhaps  even  misleading. 

7.4.  Three-Dimensionality 

The  far-wake  structure  behind  a  two-dimensional  cylinder  appears  to  be 
three-dimensional.  Our  flow- visualization  experiments  reveal  that  three- 
dimensionality  begins  around  100  diameters,  where  secondary  vortices  exhibit  a 
waviness  which  amplifies  downstream.  Waviness  of  vortex  lines  in  the  x-z  plane 
has  been  observed  previously  by  other  investigators.  Hama  (1957)  and  Tritton 
(1959)  are  two  of  the  earliest  reporters  of  this  phenomenon.  More  recently  Ger- 
rard  (1978)  and  Slaouti  and  Gerrard  (1981)  investigated  the  influence  of  flow 
non-uniformities  and  end  effects  on  vortex  lines.  Recent  numerical  work  by 
Robinson  (1984)  indicates  that  a  vortex  street  can  be  more  unstable  to  three- 
dimensional  disturbances  than  to  two-dimensional  disturbances  over  a 
significant  range  of  street  spacing  ratios.  Hama  (1963)  described  how  a  vortex 
filament  can  be  distorted  three-dim ensionaiiy  and  progressively  by  its  own 
induction.  This  mechanism,  along  with  influences  from  other  vortices  on  either 
side  of  the  street,  is  postulated  here  for  our  observed  progressive  three- 
dimensionality. 

Although  the  evidence  is  not  conclusive,  we  may  speculate  that  such  distor¬ 
tions  eventually  lead  to  "pinching  off"  of  vortices  from  either  side  of  the  wake  to 
form  loops.  As  pointed  out  by  Roshko  (1978),  loops  of  this  kind,  when  averaged 
for  a  long  time,  may  explain  Grant's  "vortex-pair  eddy"  which  has  puzzled  us  for 


so  many  years. 
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Appendix  A 


Numerical  Streak  lines 

I 

A.1.  Procedure 

Following  Hama  (1962),  Euler's  method  is  used  to  calculate  streaklines  gen- 
erated  in  a  given  velocity  field.  Particles  are  released  from  a  fixed  point  at  each 

■ 

time-step  Af .  The  position  of  each  particle  is  advanced  every  time-step,  accord¬ 
ing  to 

»  U(t+At)  =  Xit)  +  ^p-(At)2  +  ...  (A.1) 


I 


■ 


In  the  computer  code  Equation  (A.1)  is  calculated  to  second  order  for  every  par¬ 
ticle  at  each  time-step.  Connecting  a  string  of  these  particles  together  (all  of 
which  had  been  introduced  from  the  same  point)  constitutes  a  streakline. 


A.  2.  Viscously-Decaying  Vortex  Street 

The  two-dimensional  velocity  field  used  for  the  numerical  simulation  in 
Chapter  3  is  described  here.  Following  Lamb  (1945),  the  inviscid  velocity  field  of 
an  infinite  row  of  vortices,  each  with  strength  T  ,  and  coordinates  (0,0),  (±  a,0), 
(±  2a  .0)  ...  is 
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To  model  a  Karman  vortex  street  two  such  rows  are  superposed.  One  of  the  rows 
has  its  vortex  centers  midway  between  those  of  the  other  row,  with  the  circula¬ 
tion  equal  but  opposite;  the  two  rows  are  separated  by  distance  6  in  the  in¬ 
direction.  Such  a  configuration  has  a  self-induced  velocity  in  the-  negative  x- 
direction, 


Ut 


2a 


■tanh 


rr  fa 
a 


(A.  4) 


Adding  a  constant  freestream  velocity  to  the  array  of  vortices  completes  the 
in  viscid  model;  each  vortex  convects  downstream  with  velocity  Ue  =  Um  -  [/<. 


To  add  the  effects  of  viscosity,  the  vortices  are  allowed  to  diffuse  as  their  life¬ 
times  increase.  Each  shed  vortex  is  assumed  to  be  axisymmetric  about  its 
center  and  viscously  decaying,  with  its  induced  tangential  velocity  given  by 


(A.  5) 


where  r  is  the  distance  from  the  center,  T  is  the  circulation  or  strength  of  the 
vortex,  and  v  is  the  kinematic  viscosity  of  the  fluid  medium.  At  time  t  =  0 
Equation  (A. 5)  reduces  to  that  of  an  inviscid  axisymmetric  vortex.  As  time 
increases,  the  vorticity  spreads  out  (diffuses)  in  such  a  manner  that  total  circu¬ 
lation  remains  constant.  Far  enough  away,  the  induced  velocity  of  this  vortex  is 
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no  different  than  that  for  an  inviscid  rectilinear  vortex.  Thus,  to  model  viscos¬ 
ity,  the  inviscid  vortices  of  the  idealized  Karman  street  are  replaced  by  Equation 
(A. 5)  near  regions  of  interest,  i.e.,  close  to  the  expanding  core  of  the  vortex. 
Time  t  is  chosen  as  the  ’‘lifetime”  of  the  vortex;  in  the  computer  code  t  =  x/Uc 
where  x  is  the  downstream  distance  the  vortex  has  travelled  and  Uc  is  its  con¬ 
vection  speed.  Vortices  are  thus  "born”  periodically  at  x  =0,  and  are  con- 
vected  downstream  to  form  the  Karman  vortex  street  pattern. 


It  should  be  kept  in  mind  that  the  above  model  is  not  an  exact  solution  of  the 
equations  of  motion,  but  as  discussed  by  Hooker  (1936),  it  is  nevertheless  a  very 
reasonable  model.  The  constants  in  Equation  (A. 5)  were  chosen  to  simulate  the 
flow  represented  by  Figure  3.1.  In  dimensionless  variables  Equation  (A.5) 
becomes 


U-9 


~(r*)zRe 
4f  * 


(A.6) 


where 


Ug  = 


ue/um  ' 


r*  =  rM  ,  r* 


r 

Umd  ' 


Re  = 


U„d 

v 


For  the  calculations  shown  in  Figure  3.2,  T'  =  1.4,  (an  experimental  result  of  • 


Matsui  it  Okude  (1981)),  Re  =  90,  a/d.  =  6,  and  bM  =  1.69.  The  latter  two 
values  are  close  to  those  measured  in  Figure  3.1.  Bused  on  experimental  results 
of  Hooker  (1936),  each  vortex  is  given  an  “apparent  age"  at  x  =  0  of  t  *  =  2. 
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Appendix  B 

Cross-Wire  Wake  Surreys 

B.l.  Cross-Wire  Calibration  Procedure 

A  pitot-static  tube  and  electronic  manometer  were  used  to  calibrate  cross¬ 
wire  voltages  against  tunnel  velocity,  immediately  prior  to  each  run.  This  was 
done  by  simultaneously  recording  the  voltage  from  each  wire  along  with  manom¬ 
eter  output  -  the  cross-wire  probe  and  pitot-static  tube  being  positioned  near 
each  other  and  in  the  freestream  of  the  tunnel,  thirty-two  two-second  averages 
were  recorded  as  the  tunnel  velocity  was  varied.  (The  hot-wire  voltages  were 
first  conditioned  to  lie  within  the  range  0  to  10  volts,  so  as  to  use  the  full  range 
of  the  A/D  converter.)  The  voltage  from  each  wire  was  then  plotted  against  velo¬ 
city  normal  to  that  wire.  A  typical  calibration  (fitted  to  a  fourth-order  polyno¬ 
mial)  is  shown  in  Figure  B.l.  This  calibration  curve  was  subsequently  used  to 
calculate  the  instantaneous  u  and  v  velocities  for  each  sample  in  the  wake 
survey. 

B.2.  Wake  Surveys 

Because  of  the  large  (a  550)  length-to-diameter  ratio  of  the  wires,  hot-wire 
end  effects  were  neglected,  in  accord  with  the  experimental  conclusions  of 
Champagne,  Sleicher,  and  Wehrmann  (1967).  Thus,  a  simple  trigonometric 
decomposition  (the  "cosine  law")  was  used  to  obtain  u  and  v  for  each  cross- 
wire  sample;  i.e., 

U  =  (lX+45*  +■  U.  _4j»)  cos45°  ,  (B.l) 


•  i 


•  « 


•  « 


•  ( 


•  « 


•  i 


•  4 
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v  =  (u+4fi*  -  u^M.)cos45a  , 


(B.2) 


where  u+*5.  and  u^s*  are  the  velocity  components  normal  to  the  +  and  -45° 
wires,  respectively,  as  defined  in  Figure  B.2.  u  and  v  are  the  instantaneous 
velocities  which  we  now  write  in  terms  of  mean  and  fluctuating  components, 

u  =  U  +  u'  ,  (B.3) 

v  =  V  +  v'  ,  (B.4) 


•  * 


•  ( 


•  4 

For  the  case  of  circular  cylinder  wakes  at  Re  =  150  and  500,  the  specifics  of 
the  wake  surveys  are  as  follows:  The  cylinder  diameter  was  0.16  cm;  Reynolds 
numbers  of  150  and  500  were  obtained  with  freestream  velocities  of  1.52  and 
4.97  m/s  respectively.  Twelve  ar-locations  were  chosen,  namely  x/d  =  25.  50,  *  ' 

75.  100,  125,  150,  200,  250,  300,  400.  500  and  750.  At  each  of  these,  data  were 
recorded  at  thirty-two  y  -positions  (the  spacing  in  the  y  -direction  was  adjusted 
according  to  the  local  width  of  the  wake).  The  sampling  rate  was  2000  Hz  for  •  < 

each  of  the  two  wires;  the  total  number  of  samples  per  wire  at  each  position  in 
the  wake  was  8192.  Thus,  the  total  sampling  time  at  each  position  was  slightly 
more  than  four  seconds. 

The  same  cylinder  was  used  for  the  data  at  Re  =  140,  at  a  freestream  velocity 
of  1.38  m/s.  The  z-iocations  chosen  for  the  survey  were  more  compact  than 
those  above;  x/d  =  40,  50.  60,  75,  100,  125,  150,  175,  200,  250,  300,  and  350.  The  • 

porous  plate  chosen  for  detailed  cross-wire  surveys  was  a  47%  solidity  screen  of 
width  d  =  1.78  cm,  whose  wire  diameter  and  mesh  size  were  0.43  and  1.59  mm 
respectively.  Surveys  were  made  at  x/d  =  1,  2,  4,  7,  10,  15,  20,  25,  30,  40,  and  50. 
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Two  channels  of  data  were  recorded  simultaneously  -  the  voltage  output 
from  the  +45°  and  -45"  wires.  The  calibration  curve,  an  example  of  which  is 
shown  in  Figure  B.l,  was  used  to  calculate  u^s.  and  Equations  (B.l)  and 

(B.2)  provide  u  and  v,  which  we  could  average  over  the  entire  record  to 
obtain  U  and  V.  The  fluctuating  velocity  components  u'  and  t/  are  then 
derived  from  Equations  (B.3)  and  (B.4).  Figures  B.3  -  B.5  show  U  profiles  for 
various  downstream  positions  in  a  cylinder  wake  at  Re  =  140,  150,  and  500 
respectively. 

B.3.  Analysis  of  the  Data 

Let  us  consider  the  wake  survey  for  a  circular  cylinder  wake  at  Re  =  150. 
Having  recorded  u  and  v  for  each  sample,  there  are  many  ways  to  analyze  the 
data.  Firstly  the  mean  velocity  U  is  normalized  in  the  usual  manner: 

■  <b« 


where  U0  is  the  centerline  velocity,  Um  is  the  freestream  velocity,  and  y*  is 
the  normalized  transverse  coordinate, 


0 

V  = 


(B.6) 


In  Equation  (B.5),  the  constant  In 2  was  chosen  such  that  U*  -  at  y  *  =  1  or 

-1.  <5  is  the  wake  half -width;  Vo  Is  the  transverse  distance  corresponding  to  the 
centerline  of  the  wake.  U0,  6,  and  y 0  were  chosen  to  give  the  best  (least- 
squares)  fit  to  Equation  (B.5).  Figure  5.14  shows  an  example  of  how  well  the 
data  collapse  onto  this  Gaussian  wake  profile.  The  wake  half-width  6  is  an 
important  quantity  which  is  used  to  normalize  frequency  /  : 
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(B.7) 


•  4 


Another  important  quantity  is  the  centerline  velocity  defect, 

•  < 

Wa  =  l  (B.8) 

•  < 

Figures  5.15  and  5.18  provide  examples  of  how  <J  and  W0  vary  with  downstream 
distance  x/tf. 

It  is  not  certain  why  such  a  large  scatter  exists  in  these  measurements.  Our  m  ( 

conjecture  is  that  non-uniform,  unsteady  boundary  layer  transition  on  the  walls 
of  the  test  section  leads  to  a  low-frequency  oscillation  of  the  freestream. 

Ideally,  data  should  be  gathered  at  all  wake  positions  simultaneously  to  avoid  _ _ 

any  such  problems;  this  was  not  possible  however,  since  only  one  cross-wire  was 
used.  Hence  freestream  conditions  could  have  changed  during  the  course  of  a 
survey.  Four  seconds  may  not  have  been  a  long  enough  sampling  period  to  aver¬ 
age  away  these  (very-low-frequency)  fluctuations.  • 

In  addition,  the  height  (y  -direction)  of  the  cross-wire  was  approximately  the 
same  as  the  cylinder  diameter.  Hence,  close  to  the  cylinder,  where  the  wake 
width  is  small,  we  would  not  expect  good  resolution  in  measurements  of  mean  * 

velocity.  For  this  reason,  the  smoothing  in  Figures  5.15  and  5.16  is  weighted  in 
favor  of  far-wake  data. 

• 

In  all  subsequent  numerical  calculations  in  which  6  and  W0  were  required, 
smoothed  values  were  used  rather  than  the  measured  values,  as  shown  in  Fig¬ 
ures  5.15  and  5.16.  It  is  important  to  note  here  that  power  spectra,  probability 
density  functions,  and  other  quantities  based  on  the  instantaneous  fluctuations  • 
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are  not  affected  by  this  smoothing. 

There  are  several  methods  by  which  we  can  search  for  wake  structure,  given 
the  time-traces  of  u'  and  v' : 

a)  Time  plots  -  An  example  plot  of  v'  AJm  versus  sample  time  is  shown  in 
Figure  5.9.  We  have  also  plotted  W /Um,  and  the  instaneous  Reynolds  stress 
—u’v'/Um  in  Figures  B.6  and  B.7,  respectively. 

b)  Amplitude  spectrum  -  A  fast-Fourier-transform  (FFT)  algorithm  can  be 
used  to  obtain  the  amplitude  spectrum  of  the  signal.  This  gives  us  a  quanti¬ 
tative  description  of  the  fluctuation  magnitude  as  a  function  of  frequency; 
an  example  is  shown  in  Figure  5.1.  The  method  of  modified  periodograms 
was  used  to  generate  power  spectra,  as  discussed  in  Rabiner  and  Gold 
(1975).  (Amplitude  spectrum  is  defined  as  the  square  root  of  the  power 
spectrum.)  Spectra  of  the  freestream  tunnel  conditions  for  each  wake  sur¬ 
vey  are  shown  in  Figures  B.8  and  B.9. 

c)  Probability  density  -  If,  in  the  time-trace  of  the  fluctuating  velocity, 
there  exist  identifiable  "groups"  or  "bursts”  of  structure,  the  power  spec¬ 
trum  can  be  calculated  for  each  such  group  individually.  The  most  prom¬ 
inent  frequency  fp  of  each  group  is  then  determined,  and  the  results  of  the 
entire  time-trace  are  displayed  in  terms  of  probability  density  of  spectral 
peaks  as  a  function  of  frequency.  Probability  density  is  a  bin-independent 
quantity  calculated  as  follows:  frequency  is  divided  into  bins  spaced  apart 
by  A/.  We  then  count  the  number  of  times  fp  falls  into  each  bin.  The 
result  is  called  a  probability  distribution,  which  is  not  independent  of  bin- 
width  A/ .  Dividing  by  the  total  area  under  the  probability  distribution 
curve  gives  the  probability  density,  which  by  definition  has  unit  area,  and  is 
therefore  bin -independent.  An  example  of  such  a  spectral  probability 
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density  function  (PDF)  is  given  in  Figure  5.10. 

To  identify  a  group,  an  algorithm  was  written  to  scan  the  data,  searching 
for  a  place  where  peak-to-peak  fluctuations  remain  above  some  velocity 
threshold  for  a  longer  period  than  some  time  threshold.  In  the  code  the 
velocity  threshold  ut  was  chosen  to  be  vlS  u'ma  .  where  v'm,  was  calcu¬ 
lated  for  the  entire  time-trace.  The  factor  y/2  was  chosen  since  for  a  sine 
wave  v'p-j,  =  vVm*’  The  proper  time  threshold  to  choose  is  somewhat  less 
clear.  If  too  small,  the  frequency  resolution  of  the  FFT  would  be  poor;  if  too 
large,  some  of  the  shorter-lasting  groups  of  structure  would  not  be  con¬ 
sidered.  The  time  threshold  chosen  in  our  calculations  therefore  varied 
with  location  and  Reynolds  number.  Typically,  tt  was  chosen  so  as  to  allow 
at  least  two  cycles  of  a  disturbance  at  the  minimum  expected  frequency. 

The  PDF's  were  weighted  by  the  amplitude  of  fluctuations  within  a  group 
-  i.e.,  those  groups  with  large-amplitude  fluctuations  contributed  more  to 
the  PDF  than  groups  with  smaller  amplitudes. 
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Appendix  C  • 

In  viscid  linear  Stability  Calculations 

For  parallel  two-dimensional  mean  flows,  with  U=U(y)  only,  inviscid  linear  * 

hydrodynamic  stability  theory  may  be  applied.  Basically,  a  small  unsteady  per¬ 
turbation  is  superposed  on  the  known  steady  mean  solution;  analysis  of  the 
resulting  set  of  equations  can  predict  which  disturbances  are  likely  to  be  • 

amplified.  Following  Gaster  (1965),  let  the  perturbation  stream  function  be 
given  by 

=  <p(y)  (C.l)  * 


where 


“■  v  v  =  ^  <cs> 

a  is  called  the  wave  number,  while  p  is  the  frequency;  in  the  general  case  both 
a  and  p  may  be  complex.  Substituting  ^  and  the  known  mean  solution  U( y) 
into  the  equations  of  motion  yields  the  Orr-Sommerfeld  equation.  Neglecting 
viscosity,  the  resulting  equation  is  the  simpler  Rayleigh  equation, 

(C/(y)  -  p/a)(V"  -  a*<p)  =  0,  (C.3) 


where  primes  denote  differentiation  with  respect  to  y . 

For  the  present  case  of  a  far-fleld  two-dimensional  wake,  an  appropriate 
parallel  mean  solution  is  the  Gaussian:  • 
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C/(y)  =  um(i  -  r,  f). 


(C.4) 


where 


_  U(V )  -  e-OaZKv/a)* 

~  U0-U„ 


(C.5) 


(C.6) 


d  is  the  wake  half-width,  U0  is  the  wake  centerline  velocity,  and  the  constant 
ln2  is  chosen  so  that  f  at  ^-=  ±  1.  Putting  (C.4)  into  (C.3)  results  in 

(C-O(^M“<*V)  +  f>  =  0.  (C.70 


with 


C 


aUm 


(C.B) 


Kubota  (1983)  has  obtained  solutions  for  wave  speed  C,  given  a  =  real  (tem¬ 
poral  mode);  and  also  for  wave  number  a  ,  given  p  =  real  (spatial  mode).  In 
the  latter  case,  spatial  amplification  is  equal  to  -a*  ;  in  Figure  5.12  we  have  plot¬ 
ted  -a *  vs.  p  for  various  values  of  FK<j  •  For  the  temporal  mode,  the 
corresponding  spatial  growth  rate  is  approximated  by  g  =  api/pr,  according  to 
Gaster  (1965).  This  growth  rate  as  a  function  of  pr  is  plotted  in  Figure  5.13. 


-82- 


In.  Figure  C.l  the  calculations  at  =  0.3  are  compared  with  experimental 
results  of  Derm  lie  (1983).  He  measured  growth  (and  decay)  rates  of 
acoustically-excited  fluctuations  in  the  far  wakes  of  circular  cylinders  at  low 
Reynolds  numbers.  Also  shown  on  the  plot  are  recent  calculations  by  Hultgren 
(1984)  for  the  case  of  viscous  spatial  linear  stability  of  a  far  wake  at  W0  =  0.3  , 
and  at  Re4  =  70.  Re4  is  based  on  the  wake  half-width  and  velocity  defect. 


Re4  — 


V0  Um  6 
v 


(C.9) 


Reynolds  number  based  on  cylinder  diameter  d  is  also  listed  in  Figure  C.l  for 
comparison. 
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Figure  3.1:  Circular  cylinder  wake  at  Re  -  90  :  smoke-wire  at  x/d  -  a) 
b)  50.  c)  100.  and  d)  150. 


introduced  at  x/d  =  a)  0.  and  b)  100. 


gure  4.2:  Circular  cylinder  wake  at  Re  -  140 ;  smoke 
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Figure  5.1:  Amplitude  spectrum  of  v'/U„  at  inflection  points  of  mean 
velocity  profile;  circular  cylinder  wake  at  Re  *=  150,  x/d'=25. 
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Figure  5.3:  Amplitude  spectra  at  several  downstream  distances; 
circular  cylinder  wake  at  Re  -150. 


cylinder  wake  at  Re -150. 


Figure  5.5:  Amplitude  spectra  at  several  downstream  distances; 
circular  cylinder  wake  at  Re  -140. 


cylinder  wake  at  Re -  440. 


Figure  5.7a)  :  Amplitude  spectra  at  several  downstream  distances; 
circular  cylinder  wake  at  Be  -500. 


Figure  5.7b)  :  Normalized  amplitude  spectra  at  several  downstream 
distances;  circular  cylinder  wake  at  Re  -500. 
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Figure  5.7c)  :  Fully  normalized  amplitude  spectra  at  several  downstream 
distances;  circular  cylinder  wake  at  Re  -500. 


Figure  5.8a)  :  Growth-decay  cycles  of  selected  discrete  frequencies; 
circular  cylinder  wake  at  Re  -500. 


Figure  5.8b)  :  Normalized  growth-decay  cycles  of  selected  discrete 
frequencies;  circular  cylinder  wake  at  Re -500. 


Figure  5.10:  Probability  density  functions  of  v'  fluctuations; 
circular  cylinder  wake  at  Re  -150. 


Figure  5.11:  Average  number  of  cycles  per  group  versus  downstream 
distance;  circular  cylinder  wake  at  Re -150. 
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Figure  5.12:  Spatial  inviscid  stability  solutions.  Gaussian  wake  profile. 


Figure  5.13:  Temporal  inviscid  stability  solutions.  Gaussian  wake  profile. 


Figure  5.14:  Mean  velocity  profiles:  circular  cylinder  wake  at  Re  «  150. 
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Figure  5.15:  Wake  half-width  thickness  vs.  downstream  distance: 
circular  cylinder  wake  at  Re  -  150. 
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Figure  5. IB:  Velocity  defect  vs.  downstream  distance; 

circular  cylinder  wake  at  Re-  150. 
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Figure  5.17:  Numerical  simulation  of  growth-decay  cycles: 
circular  cylinder  wake  at  Re-  150. 
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Figure  5.18:  Peak  frequency  vs.  downstream  distance:  circular 
cylinder  wake  at  Re  -  150. 
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Figure  5.19:  Peak  frequency  vs.  downstream  distance;  circular 
cylinder  wake  at  Re  -140. 
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Figure  5.20:  Peak  frequency  vs.  downstream  distance;  circular 
cylinder  wake  at  Re  -500. 
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Figure  5.21:  Amplitude  spectra  at  several 
y-locations  :  cylinder  wake  at 
Re  -  300,  and  at  x/d-750. 
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Figure  5.22:  Prominent  frequency  versus  y;  circular  cylinder  wake  at 
Re  -  300,  and  at  x/d  -  750 . 


a)  high  frequency,  b)  low  frequency  structure 
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Figure  5.24:  Simultaneous  smoke-wire  and  cross-wire  records,  circular 
cylinder  wake  at  Re  -  150 ;  smoke-wire  at  x/d  -  275,  cross¬ 
wire  at  x/d  -  400,  y/d--4. 


-107- 


N 

X 


>. 

o 

c 

M 

3 

O 

© 

C- 

u. 


in  o  m  o  in  o 

C\J  (U  "rt  Ti 

apn^Ttdujv  tej^oads  "n/,A 


* 


» 


» 


» 


Figure  5.26:  Amplitude  spectra  at  several  downstream  distances: 

porous  flat  plate  wake  (o  -  47%)  at  Re  -5000. 
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plate  wake  (a  -  47%)  at  Re  -5000. 


re  at  x/d  =■  a)  a,  b)  100,  and  c)  200,  and  at  y/, 


...  •  r-kcrk • 

f*. -**•'•'?>  ,, 

^<5>  if<^-  '(r*'u'WfnrA»i 

rjiv  f;r:^j 

'^t  , 

•••  '**-"rir;  !?((&(  ^inuPKnif 

■('■pry-  t  lj«t 
*r  I’’*'  ’*  f*"5^  */3!I* 

;  K'i  it  •»•*♦;/;*' ' .: x v<t,tt r.tHUmtnffMih 

rf.?  <1  iwjfc  *«  ( J& 

,  vi » ’.ffj?yh <a4 

v**> utu^j «« 

<3/^(,«Ar*  if  "• 

^  »** 

.  w,„s  ♦  *«•»  t»> »•* ►*  >1  t"  tlll'lt  Hljii  jj|f 

o»  ,.,AV»  "i  i*i|»  u*  Ml11  *"ilMl!||! 

\  r  t<*  )  ill i  o  ii|  I  j  ill  il;i 


i 

i 


Figure  6.4:  Edge  and  plan  views  of  a  circular  cylinder  wake  at  Re  =  190 ; 

smoke-wire  at  x/d  =  a)  8,  b)  100,  and  c)  200,  and  at  y/d  =  0. 


Figure  B.l:  Cross-wire  calibration,  polynomial  fit  to  order 
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Figure  B.3:  Mean  velocity  profiles  for  circular 
cylinder  wake  at  Re -140. 
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Figure  B.5:  Mean  velocity  profiles  for  circular 
cylinder  wake  at  Re -500. 
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gure  B.6:  Time-trace  of  u'/U„;  circular  cylinder  wake 
at  Re  -  150,  x/d  -  200,  and  y/d“1.5. 


Figure  B.7:  Time-trace  of  Reynolds  stress,  -u'vVlU2;  circular 
cylinder  wake  at  Re  -  150,  x/d  -  200,  and  y/d  =  1.5. 
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Figure  B.B:  Freestream  spectra  at  tunnel  conditions  corresponding 
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Figure  B.9:  Freestream  spectra  at  tunnel  conditions  corresponding 


CALCULATIONS  AT  Wo-0.3 
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Figure  C.l:  Comparison  of  far-wake  stability  calculations  with 
experiment;  Wn*0.3. 


